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Abstract 


Lie group methods are used to find both exact and numerical similarity solutions 
for compressible perturbations to an incompressible, two-dimensional, axisymmetric 
vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for 
which the solutions are a set of two-dimensional, self-similar, incompressible vortices. 
These solutions are augmented by deriving a conserved quantity for each eigenvalue, 
and identifying a Lie group which leaves the reference flow equations invariant. The 
partial differential equations governing the compressible perturbations to these ref- 
erence flows are also invariant under the action of the same group. The similarity 
variables found with this group are used to determine the decay rates of the velocities 
and thermodynamic variables in the self-similar flows, and to reduce the governing 
partial differential equations to a set of ordinary differential equations. The ODE’s 
are solved analytically and numerically for a Taylor vortex reference flow, and nu- 
merically for an Oseen vortex reference flow. The solutions are used to examine the 
dependencies of the temperature, density, entropy, dissipation and radial velocity on 
the Prandtl number. Also, experimental data on compressible free vortex flow are 
compared to the analytical results, the evolution of vortices from initial states which 
are not self-similar is discussed, and the energy transfer in a slightly-compressible 
vortex is considered. 
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Chapter 1 
Introduction 


A compressible flow is one in which a change in pressure over some characteristic 
length scale of the flow results in a corresponding change in density. Similarly, a com- 
pressible vortex is a vortex containing large tangential velocities which create strong 
pressure gradients; these pressure gradients, in turn, produce substantial density vari- 
ations across the vortex. The main objective of this thesis is to present analytical self- 
similar solutions of the equations describing the viscous decay of a two-dimensional, 
axisymmetric, slightly-compressible, free vortex. A quantitative characterization of 
the problem at hand is given in § 1.1. However, let us first consider the current state 
of basic research on this subject. 

Vortices have long fascinated fluid dynamicists and although the total number of 
investigations involving vortices is quite large, there has been relatively little work 
on compressible vortex flow. For example, a search of the citations in Stanford Uni- 
versity’s science journal holdings reveals that, since 1969, the university has collected 
about 18, 500 articles on vortices. A similar search shows that, during the same 
time period, the university libraries acquired around 330 compressible vortex articles. 
Assuming these numbers are representative of the general state of vortex research, 
slightly less than 2% of the work on vortex flow takes into account the effects of com- 
pressibility. Further, studies of flows containing compressible vortices do not usually 
focus on the basic physics or structure of a vortex, but concentrate on other aspects 
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of the flow such as shock/vortex interaction, blade/vortex interaction, or vortex 
breakdown. 

The dearth of basic research on compressible vortices is especially surprising given 
the fact that they play a key role in many fluid flows of physical interest: turbulent 
combustion in piston and turbine engines, strong wing-tip vortices, flows with aerody- 
namically generated noise, shock diffraction around sharp corners, and astrophysical 
flows such as those occurring in the solar corona or in the gaseous disks of spiral 
galaxies. Most fluid dynamicists are surprised to find that there are no simple an- 
alytical solutions describing the viscous decay of a compressible free vortex. The 
mathematical solutions to the equations of compressible vortex flow are complicated 
and must usually be numerically obtained. The primary goal of this thesis is to de- 
scribe a set of simple analytical solutions for compressible free vortex flow. This work 
is motivated by the fact that such solutions can yield a great amount of insight into 
the fundamental physical behavior of compressible vortices, as well as provide a quick 
and easy means of validating computational and experimental results. 

The next section further describes the flow under consideration and is followed 
by a brief literature survey of work on compressible vortices and an overview of the 
thesis. 


1.1 Orientation 

The lack of work on compressible free vortices is partly due to their complexity. In 
contrast, the two-dimensional motion of an isolated free vortex in an incompressible 
flow is well understood. Conservation of mass requires that there is no radial velocity 
in the vortex, and the Reynolds number Re, governs the vortex’s evolution. For very 
large Reynolds numbers, Re oo, the flow is essentially inviscid, and the vortex 
motion is nearly steady. Otherwise, for Re ~ 0(1) or smaller, one finds that the 
viscous diffusion of angular momentum is the process which dominates the flow. 

When the compressibility of fluid within the vortex is important, the flow has 
an analogous dependence on the Reynolds number, but analyzing the flow becomes 
more complicated. In studying a two-dimensional, compressible vortex one must 
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examine the effects of heat conduction, viscous dissipation, compressibility, and radial 
convection on the motion and structure of the vortex [7]. In addition to Re, solutions 
for compressible vortex flow require that four other parameters be considered: the 
Mach number M , the Prandtl number Pr , and the temperature dependent values of 
viscosity fi, and thermal conductivity k,. 

To get a quantitative feel for the nature of the compressible vortices under investi- 
gation, consider the vortex pictured in figure l.l(a-e). In general, the radial profiles 
of the velocities and thermodynamic variables (pressure, density, temperature, etc.) 
in such a vortex may be slightly different, but the sketch illustrates the main struc- 
tural features of a compressible vortex. Here r is the radial coordinate and the origin 
of the system is the center of the vortex (see figure 1.2). Since the present work is 
limited to two-dimensional axisymmetric vortices, the flow does not vary in the tan- 
gential or axial directions. Here we will also assume that the far-field (r oo) flow 
conditions are constant. 

The maximum tangential velocity v m occurs at the radial location k (figure 1.1a); 
k will be called the core size of the vortex. In the following analysis, k and v m are the 
length scale and the velocity scale, respectively. Together with the far-field viscosity, 
density, and sound speed, k and v m are used to form the Reynolds number and Mach 
number of the vortex. 

As a rule of thumb, the maximum tangential velocity can be considered u large” 
when its magnitude is greater than or equal to about 30% of the far-field speed of 
sound. Thus, for a Mach number of M ~ 0.30 the vortex will be considered com- 
pressible; when the Mach number is lower than M & 0.3 a simplified analysis which 
assumes the flow is completely incompressible should generally suffice. However, the 
vortex studied here is slightly- compressible and so v m cannot be too large either. As 
will be further explained in §2.2 the present study uses the Rayleigh- Janzen expansion 
[38] to derive equations describing compressible perturbations to an incompressible 
vortex reference flow. Colonius, Lele, & Moin [12] compare the results of computations 
of the full Navier-Stokes equations and numerical solutions for a slightly compress- 
ible vortex flow (found using the Rayleigh- Janzen expansion). They show that the 
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approximation is valid for Mach numbers as large as M « 0.67. The same flow con- 
ditions apply to the vortices in this work, so we will take the range of Mach numbers 
under consideration to be 0.30 ~ M ~ 0.67. 

The centrifugal force in the rotating vortex flow produces a pressure minimum at 
the center of the vortex. The corresponding density variation in the vortex is regarded 
as significant when the non-dimensional depth of the density well A p*, which is the 
difference between the density in the center of the vortex and the density in the 
far-field p^, normalized by the far-field density, is A p* ~ 0.3. 

Consider the vortex studied by Mandella & Bershader [26]. A two-dimensional, 
compressible, free vortex was created in a shock tube by the diffraction of a shock 
around a 3.7 cm chord-length NACA 0018 airfoil (see §3.3). The strength of the in- 
cident shock was M s « 1.8, and the airfoil’s angle of attack was 30°. The vortex was 
embedded in a free stream with a Mach number of M , « 0.52 and a Reynolds num- 
ber, based on airfoil chord-length, of Re^ « 4.4 x 10 5 . Holographic interferometry 
and pressure transducers were used to measure the density and pressure, respec- 
tively, in the vortex studied. The measured density and pressure distributions in 
the vortex and the estimated tangential velocity in the flow give k « 9.3 x 10 -1 mm, 
v m « 2.3 x 10 2 m s" 1 , M a 6.7 x 10 -1 , and A p* « 5.5 x 10” 1 . The vortex Reynolds 
number (based on v m , and a far-field kinematic viscosity of « 1.1 x 10~ 5 ) is 
Re ze 2.0 x 10 4 . The magnitude of the radial velocity in a compressible free vortex 
is generally very small (figure l.le). Mandella & Bershader estimate the magnitude of 
the maximum radial velocity in the vortex they studied to be (max |u|)/u m « 1.3 x 10“ 3 
However, Colonius et al. [12] have shown that if the far-field flow conditions are as- 
sumed to be constant, Mandella’s estimate of the radial velocity is too large and 
therefore violates the Second Law of Thermodynamics. They suggest the radial ve- 
locity should be around two orders of magnitude smaller. Nevertheless, the vortex in 
Mandella & Bershader’s experiments clearly exhibits the tangential velocity, density, 
and pressure profiles characteristic of a compressible vortex. The issues of the entropy 
distribution and the entropy production in a compressible vortex will be discussed in 
§§ 6.3 and 6.4. 








CHAPTER 1. INTRODUCTION 


7 



Figure 1.1: The structure of a two-dimensional, compressible, free vortex — radial 
profiles of (a) tangential velocity, (6) pressure, (c) density, (d) temperature, (e) radial 
velocity. 

The vortices considered in this work are called free vortices because there are 
no boundaries in the flow for 0 < rfk < oo. Therefore, the fluid’s motion is only 
determined by its internal shear stresses and the exhange of momentum between 
neighboring fluid particles. This is in contrast to a driven vortex in which the move- 
ment is influenced by an externally controlled surface — such as a rotating cylinder 
centered at the origin of the flow [6], [23]. Additionally, we will assume that there are 
no acoustic waves present in the flow (see §2.2.1). 

Note that in this work the interferograms and data obtained in the experiments 
of Mandella [25] will be used to describe the formation of compressible free vortices, 
study their entropy distributions, and to construct initial conditions for our self- 
similar vortex solutions. However, one must be careful to remember that the vortex 
solutions presented in this theoretical study correspond to idealized flow conditions. 
The flow-field in the experiments is less “hospitable” and less well-characterized. The 
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Figure 1.2: Sketch of the coordinate system and flow-field. Circular lines depict 
the rotating fluid in the vortex. The curved line radially emanating from the origin 
schematically represents the tangential velocity profile of a Taylor vortex [36]. 


timescale of diffusion for the vortex studied by Mandella would be t d ss 7.9 x 10~ 2 s, 
based on a coresize of k « 9.3 x lO -1 mm, a viscosity of p « 1.9 x 10 -5 Nsm' 2 , and 
a density of ~ 1.8 kgm 3 ; the observation period in Mandella’s experiments 1 , 
tabs ~ 500 ps, is about two orders of magnitude smaller than t d . In comparison, the 
acoustic timescale for the vortex is t a & 2.7 ps based on the same value of k and a far- 
field sound speed of a 0 0 « 3.4 x 10 2 ms _1 . At early times in the motion of the vortex 
in Mandella’s experiment, the presence of acoustic waves in the flow may affect the 
vortex’s evolution. Figures 1.3 and 1.4 show a compressible free vortex which has been 
created using the same conditions as the vortices studied by Mandella & Bershader 
[26]. Notice the weak circular shockwave, the non-linear acoustic waves, and the rolled 
up vortex sheet surrounding the vortex. The experimental flow is quite complex and 
the far-field flow conditions may be time dependent 2 . The acoustic waves propagate 
away from the vortex just after the vortex has formed. In this experiment, and in 

1 During this observation period, the free vortex travels a distance of « 9.0 cm and executes about 
130 rotations or “eddy turnovers” . 

2 Over the short observation time of Mandella’s experiment, the distributions of radial veloc- 

ity and entropy may be very sensitive to small changes in the far-field pressure and density; see 
von Ellenrieder, Kao, &; Bershader [40] and Appendix A. 
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the experiments of Mandella & Bershader the observations are made within 500 ^secs 
of the vortex’s formation. Because of the presence of these waves, and the short 
observation period of Mandella’s experiments, a direct comparison cannot be made 
between the time evolution of the vortex studied by Mandella and the viscous decay 
of the self-similar vortex solutions developed here. 



Figure 1.3: A shadowgraph of the flow-field around a vortex created under the same 
conditions used in the experiments of [26]. There are several acoustic waves present 
in the flow (circular waves centered near the vortex center) . Also, one can see several 
contact surfaces and a slip-line, which extends from the trailing-edge of the airfoil, 
being wrapped around the vortex. 


The final approximations made in the following analysis are that the coefficients 
of viscosity and thermal conductivity are constant. In summary, the following as- 
sumptions have been made: The vortex is a two-dimensional, axisymmetric, slightly- 
compressible, viscously-decaying, free vortex; the far-field flow conditions and the 
coefficients of viscosity and thermal conductivity are constant. 
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Figure 1.4: Two-color (yellow & blue) schlieren image of the flow-field around a 
vortex created under the same conditions used in the experiments of [26]. This image 
corresponds to a slightly later experimental time than shown in figure 1.3. The density 
variations in the separation bubble and vortex at the nose of the airfoil, as well as 
the shock system, are more apparent in this image than in the preceding picture. 
However, the shadowgraph in figure 1.3 more clearly captures the slipstreams present 
in the flow. 


1.2 Literature survey 

This section discusses only previous work on viscous two-dimensional free vortices. 
The reader interested in steady inviscid vortices [11], [2], driven vortex flow [23], 
[6], three-dimensional compressible vortices [27], or in compressible leading-edge vor- 
tices [9] is referred to the cited works and to the fairly comprehensive review written 
by Bershader [7]. 

Mandella [25] and Mandella & Bershader [26] discuss a compressible free vortex 
which is produced by shock diffraction in a shock tube. Interferometric density mea- 
surements are used to develop empirical formulae for the evolution of the density 
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and radial velocity in the vortex. The pressure in this vortex is measured at one 
location along its trajectory. Using these data, and assuming that the radial pressure 
gradient is balanced by centrifugal forces, the distributions of temperature and tan- 
gential velocity in the vortex are calculated. Mandella & Bershader conclude that: 
the compressible free vortex studied in their experiments is not homentropic; the ra- 
dial velocities present in the vortex arise to counteract the viscous diffusion of angular 
momentum 3 ; and the primary effects of compressibility are the inward convection and 
compression of rotating fluid. 

Sibulkin [33] investigates the transfer of total energy within a decaying, two- 
dimensional, viscous rotating flow. Low Mach number flows of both liquids and 
gases are considered. Analytical solutions for the perturbations to the total energy 
in a Taylor vortex [36] reference flow are given for early times in the evolution of 
the vortex (when Pr — 1 and 7 = 2). Sibulkin also considers the Pr variation of 
the driving terms in the perturbed energy equation. Unfortunately, in Sibulkin’s 
approach, velocity perturbations which are on the order of the Mach number squared 
are neglected in the analysis, and so the radial velocity in the vortex is assumed to be 
zero. This leads to a contradiction in the overall solution. The solution predicts that 
the total mass in the vortex increases in time, yet aside from radial convection there is 
no physical mechanism to cause the accretion of mass in the vortex. Additionally, with 
the omission of the radial and tangential velocity perturbations Sibulkin’s solutions 
for the total energy are incomplete. 

Colonius et al. [12] perform an analytical and numerical study of compressible free 
vortices. Viscosity and heat conductivity are taken into account with the simplifying 
assumption that both are constant. The Navier-Stokes equations are expanded in 
powers of Re~ l to determine the evolution of the flow on a viscous and on an acoustic 
time-scale. An ordinary differential equation (ODE) for the vortex’s radial velocity 
is found using the viscous timescale equations. Estimates of the radial velocity for 
the vortex described in the experiments of Mandella & Bershader are made using the 
ODE. The estimates are compared to both the empirical formulas mentioned above 

- - $ 

3 See appendix A for a further discussion of angular momentum conservation in a compressible 
vortex. 
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and to computational solutions of the full Navier-Stokes equations. Colonius et al. 
find that the experimentally determined radial velocity is too large, and that the 
empirical formulas violate the Second Law of Thermodynamics — the net amount of 
entropy generated in the flow would be negative. Additionally, the viscous timescale 
equations are expanded in powers of Mach number squared (M 2 ) to derive a set of 
equations describing compressible perturbations to an incompressible reference flow. 
Solutions for the case of an Oseen vortex reference flow are determined using numerical 
integration. Two different initial conditions, which are not self-similar, are used for 
the calculations. In the first case the vortex initially has a constant density; in the 
second case the vortex is initially homentropic. Colonius et al. find that the radial 
velocity in a slightly-compressible vortex is generated by viscous effects and that the 
vortex is compressed, regardless of the initial distributions of density and entropy. 

1.3 Overview 

This thesis presents a set of self-similar solutions for the compressible perturbations 
to an incompressible, two-dimensional, axisymmetric vortex. A discussion of the per- 
turbation equations is given chapter 2 and is followed in chapter 3 by a description 
of the reference flow solutions, the corresponding conserved quantities, and a dis- 
cussion of how compressible free vortices may be formed. Similarity variables which 
transform both the reference flow equations and the perturbation equations from par- 
tial differential equations (PDE’s) to ODE’s are derived in chapter 4, and some of 
the properties of the new equation set are discussed. In chapter 5 exact solutions 
and numerical solutions for the Taylor vortex [36] temperature, density and entropy 
perturbations are presented. Numerical solutions for the Oseen vortex [30] and the 
Taylor vortex base flows are compared and the influence of the Prandtl number on the 
temperature, density, entropy, and radial velocity distributions is also investigated. 
In chapter 6 we discuss the dissipation produced in the Taylor vortex, examine the 
decay of vortices from initial conditions which are not self-similar, and investigate the 
nature of the radial velocity and the energy transfer in a slightly compressible vortex. 
A brief summary of this thesis and some concluding remarks are given in chapter 7. 



Chapter 2 

Similarity and the Navier-Stokes 
Equations 


2.1 Two-dimensional, axisymmetric, 
Navier-Stokes equations 


For a two-dimensional, cylindrically-symmetrical fluid flow, the conservation equa- 
tions for mass, radial momentum, tangential momentum, and internal energy, respec- 
tively, are 

dp Id, 

m + rfr {pur) = 0 ’ ( 2 - la ) 


du du 1 dp v 

dt dr + p dr r 


2 2p d 
+ — 


1 d 

,-—{ru) 
p dr [r dr 


+ 


d_ 

dr 


2p\ld 
* ~ Tj r¥r {rU) 


dv dv uv p d 

dt U dr r pdr 


Id 

rfr {rV) 


de de p d . , „ 1 d ( dT' 

P— + PUfr = T— (™) + T + “ I Kr " 


dt 


r dr 


r dr 


dr 


(2-1 b) 
(2.1c) 
(2. Id) 


The radial coordinate r, is the distance from the center of the vortex, and t represents 
time. The radial and tangential velocities u and v (respectively), as well as the 
thermodynamic variables (internal energy e, pressure p, density p, and temperature T) 
are functions of r and t. The bulk viscosity the shear viscosity p, the thermal 
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conductivity k , and the ratio of specific heats 7, are assumed constant in this analysis. 
Additionally, the fluid under consideration is taken to be an ideal gas with equation 
of state p = pRT , where R is the specific (molar) gas constant. The term T in the 
energy equation is the dissipation function given by 


T = 2p 



/rA 2 r 2 d 
VrJ + 2 dr 




(2.2) 


The quantity T is positive definite and corresponds to the part of the viscous work 
used to deform (not accelerate) a fluid particle [37]. 

The fluid under consideration is thermally and calorically perfect, so the internal 
energy is linearly proportional to the temperature: e = C V T. C v is the the specific 
heat at constant volume, and is related to the specific heat at constant pressure C p 
by R = C p - C v and alternatively by 7 = C p /C v . Using these relations, (2.1a), and 
the equation of state, the internal energy equation (2. Id) can be written in terms of 
the pressure 


dp 

dt 


+ .| +? | ( „ )= ( 7 - I , rt fci)|(^). ( , 3) 


2.1.1 Simplifications 

The large number of parameters and the nonlinearity of the equations shown above 
makes them difficult to solve analytically. In general, analytical solutions require sub- 
stantial simplification of the boundary conditions, the initial conditions, and the form 
of the equations. For example, Chiocchia [11] and Ardalan et al. [2] assumed steady, 
inviscid flow and applied hodograph transformations to solve ideal-fluid, compress- 
ible vortex problems. Earlier, Mack [23] and Bellamy- Knights [6] investigated steady, 
viscous, compressible vortex flows. In both of the latter two works, the boundary con- 
ditions at the center of the vortex were simplified by driving the flow with a rotating 
circular cylinder. 
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In order to describe unsteady, viscous, compressible vortex flow, Colonius et al. [12] 
made several simplifications to (2.1a — c) and (2.3). The initial conditions were spec- 
ified at all radial locations by assuming that the flow starts in a state where the 
vortex has either a constant density distribution or a constant entropy distribution. 
The complexity of the equations is reduced by assuming a high Reynolds number flow 
and using low Mach number perturbation expansions. 

Among techniques for finding analytical solutions to compressible vortex flow, the 
search for symmetry under a Lie group seems like a natural approach 1 . But note that 
for a free vortex in an unbounded domain the far-field pressure, density, and tem- 
perature asymptotically converge to non-zero constant values. However, similarity 
solutions require the system of governing equations to be invariant under transla- 
tions of the thermodynamic variables. For example consider the following boundary 
condition for the far-field density: 


Um /» = A». 


where Poo is a constant. A simple test to determine if a given transformation leaves 
an equation invariant is to re-write the equation in terms of the transformed variable. 
If the equation reads exactly the same in the new variables, then the equation is 
invariant under the transformation. The boundary condition above is invariant under 
the following translation group 

P = P ~ Pool ( 2 - 4 ) 

if p' can be found such that 

lim p' — 0. 

r— too 

But, on substituting the transformation (2.4) into (2.1a), we see that the transfor- 
mation does not leave the continuity equation invariant. An extra term containing 
Poo appears in the transformed equation, but is not present in the original equation 


dp ' Id,,. 1 d . . 

~kt + ~tt\P ur ) + --z-(poour) = 0 . 
dt r dr v 7 r dr K 7 


1 A brief introduction to Lie groups is given in appendix B. 
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Since the full Navier-Stokes equations are not invariant under translations of the 
thermodynamic variables they will not admit similarity solutions for compressible free 
vortex flow; this is the main reason why analytical solutions for fully-compressible 
free vortices do not exist. However, as shown in chapters 3 and 4, the perturbation 
equations and the incompressible reference flows presented by Colonius et al. [12] are 
invariant under a three-parameter Lie group. This group allows a rich set of similarity 
solutions to be found for a slightly compressible vortex. In general, similarity trans- 
formations are used to decrease *the number of dimensions in a problem by combining 
independent variables [8], [10]. This provides two main benefits: 1) The number of 
initial or boundary conditions required to solve the problem is reduced. For the vor- 
tex under consideration, the use of similarity transformations makes it unnecessary 
to explicitly assign initial conditions for each flow variable at every point in the flow. 
For example, if the temperature at the center of a self-similar vortex is specified at 
a given time, the value of the temperature at all other radial locations (for the same 
time) is fixed. 2) The PDE’s describing compressible perturbations of the vortex are 
reduced to ODE’s which are, in this case, relatively easy to solve. 

2.2 Non-dimensional forms and 
perturbation expansions 

The similarity solutions for a slightly compressible vortex are developed from the per- 
turbation equations formulated by Colonius et al. [12]. We summarize the derivation 
of these equations for completeness, and make some additional remarks concerning 
the bulk viscosity, initial conditions, and boundary conditions. 

2.2.1 The viscous timescale equations 

Solutions of (2.1 a-c) and (2.3) vary on three time scales: viscous, convective, and 
acoustic. If it is assumed that the radial velocity in the vortex arises to counteract the 
viscous diffusion of angular momentum, the convective and the diffusive timescales 
are equivalent. The flow variables are written as functions of an acoustic timescale 
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£, a slower viscous timescale r, where dr = d t/Re, and a non-dimensional radial 
coordinate r. The flow quantities in (2.1a-c) and (2.3) are then expanded in terms 
of the viscous time r, and the resulting expressions are collected in powers of Re” 1 . 
The non-dimensional parameters used in the timescale expansions are 


T" 2 a oot ^ U 

r - t= -T-i v = — 

H H ^oo 


u „ p p 

U ~ —i P=~ , P= 2 

a oo Poo Poo® or 


rp CpT Poo ®CQ h ilC p 


at 


P 


K 


(2.5) 


h is the core size of the vortex’s initial velocity profile. Here, we define the core size 
as the radial position of maximum tangential velocity. The specific heat at constant 
pressure is C p , the far-field speed of sound is and all far-field flow conditions are 
denoted by the subscript oo. 

Under the further assumption that u ~ 0(1 /Re) and with u* = uRe , the viscous 
timescale equations are found to be 


dp Id 

Tr + -fdf {pUr) = ^ 


dp _ pv 2 
df f ’ 

A dv pu* d _ d 
P Yr + —Tf {rv) = df 


dp,.*dp 7P d 

!h + df + Tdf^ ) = 

The equation of state has the form: 


; d fv 
dr If 


1 d 

fdf {rv) 


+ ( 7 -n i g ( t ?L 

Pr f df \ df 


(2.6 a) 

( 2 . 66 ) 

(2.6c) 

(2.6<f) 


(7 — 1) - 

p = — -pT. 


(2.6e) 


7 
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Figure 2.1: Acoustic and viscous (shaded area) flow regions; a^At » ( vAt ) lf 2 . 


Assuming that there are no acoustic waves propagating into the region of interest, 
the Re~ l timescale expansion has the effect of separating the events occurring on 
an acoustic timescale from those events which occur over a longer viscous timescale. 
To see this, compare the acoustic propagation of a disturbance starting at r — 0 
and the viscous diffusion of a disturbance also starting at the origin. During a time 
At the acoustic disturbance travels a distance « a^At, and the viscously diffusing 
disturbance travels a distance « (z/At) 1/2 , where v — fi/poo. The velocity of the 
acoustic wave is around a oo, whereas the time dependent velocity of the diffusing 
disturbance is about (vj2A i) 1 / 2 . If we assume that at very early times At < v/a 2 ^ 
the fluid can be treated as a continuum 2 , the diffusing disturbance travels faster than 


2 Strictly speaking, one cannot treat the fluid as a continuum at such short times. For air at 
standard atmospheric temperature and pressure vja ^ « 1.3 x 10 -10 s. In comparison, the time 
between successive collisions of the air molecules would be around 1.4 x 10 -10 s. So for At < vja^ 
the gas will not reach thermal equilibrium after the disturbance’s occurence. To properly model the 
situation, one would have to invoke the kinetic theory of gases; for example see the work concerning 
an impulsively-started plate performed by Howarth [17] and Yang h Lees [42]. 
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the acoustic disturbance. However, here we are only interested in times At v/a ^ 

when < 2 oo At » ( vAt ) 1 ! 2 and the acoustic wave is traveling much more rapidly than 
the diffusing disturbance. For these times the viscous region of the flow is confined 
to a relatively small area which contains no acoustic waves 3 (see figure 2.1). Using 
k as a lengthscale, the dimensionless viscous time is given by At = vAt/l\ and 
the dimensionless acoustic time is given by At = a^At/k. The ratio of these two 
timescales is 

At _ vAt/lf v 1 

At a^At/li a^li Re 

This result shows that, for Re 1, an event occurring over a time At seems much 
longer to an observer measuring time on the acoustic timescale than it does to an 
observer measuring events based on a viscous timescale. For this reason the viscous 
timescale equations (2.6 a-e) are stationary with respect to the acoustic time t (there 
are no d t derivatives) and the relevant velocity scale for Re is a 

Another approximation made in the derivation of (2.66, d) is that we have assumed 
that the terms involving (/j„ -2p/3) in (2.16) and (2.3) are 0(1/ Re 2 ) or smaller. For 
dilute monatomic gases (ii u /p fa 0) and air (p,„/p, « 0.6) [37] this assumption is 
valid. However, for gases such as C0 2 ([iv/p « 1000) and H 2 (li v /p « 32) one must 
be careful in applying the approximation. 

Using the viscous timescale equations all of the necessary boundary conditions 
can be specified to 0(1/ Re). Since the flow is symmetrical about the origin, the 
tangential velocity, radial velocity, and heat flux vanish at the center of the vortex. 
Therefore, as can be shown with (2.66), the radial pressure gradient must be zero at 
the origin. Since both the temperature and pressure gradients are zero at the center 
of the vortex, the ideal gas law requires the density gradient to be zero at f — 0 as 
well. The flow domain is unbounded, and we will assume that both components of 
velocity vanish in the far-field. The pressure, temperature, and density are assumed 
to approach the constant values poo, T*,, and p <*,, respectively, as f — > oo. 

3 Note that Lagerstrom [22], Stewartson [34], Howarth [17], and Van Dyke [39] have investigated 
flows in which both acoustic waves and heat/ viscous diffusion are simultaneously present. 
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2.2.2 The low Mach number approximation 

To derive the equations for the compressible perturbations to an incompressible 
vortex, the viscous timescale equations are first re-written using the following non- 
dimensional parameters 


~ ^ 7 Vm! ~ V 

r = P < = T’ V = T 


U ~ P P ~ Poo 

u = —> P = P = ^~ J T 

v m Poo PooVm 


T = 


^ _ PcoVmk _ pC p 


P 


K 


M = 


_ t_ _ J_ 

Re Re 


(2.7) 


v m is a reference velocity — the maximum tangential velocity in the vortex’s initial 
velocity profile. The term p is similar to the standard definition of the dynamic 
pressure coefficient in aerodynamics [19] and represents the normalized deviation of 
the local pressure p from the far-field pressure p^. Note that the transformation from 
P to p is similar to the example transformation from p to p' shown in (2.4) above, and 
is simply a translation in pressure. 

The equations which are obtained after applying the change of variables in (2.7) 
are expanded in powers of M 2 with the Rayleigh- Janzen expansion 4 . In this method 
the dependent flow variables are approximated as 


/ — /o + M 2 fi + 0(M 4 ), 


( 2 . 8 ) 


4 This expansion has been used by many researchers for examining low Mach number compressible 
flow (Van Dyke [381, Lagerstrom [22]). Several combustion studies [28], [24], use a slightly modified 
form of (2.8) in which the expansion parameter is 7 M 2 . A justification for expanding the flow 
variables in terms of M 2 , rather than in odd powers of M, is given by Mahalingam [24]. He shows 
that asM-4 0, M 2 appears as a small parameter in the equations for the relative changes in velocity, 
pressure, density, and temperature in a one dimensional inviscid flow. 
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where / may represent any of the dependent flow variables. Terms of 0(1) are 
collected to yield equations describing an incompressible reference flow, and 0(M 2 ) 
terms are collected to give expressions for the compressible perturbations to the refer- 
ence flow. The reference flow terms are designated with a subscript 0 and the 0(M 2 ) 
perturbations are denoted by the subscript 1. 

In the reference flow the radial velocity is necessarily zero up = 0, the density and 
temperature are uniform p 0 = T 0 = 1, the pressure gradient balances the centrifugal 
forces created by the tangential velocity, and the tangential velocity is governed by a 
diffusion equation 


dpp _ Vp 
dr r 


(2.9 a) 


dv 0 

dr 


d_ 

dr 


l§z(rv 0 ) 


The corresponding boundary conditions are 


(2.96) 


f — 0 




(2.10a) 


and 


r — » co 


f 

Vo 0, 

< 

Po o. 


(2.106) 


To 0(M 2 ) the equations for conservation of mass, radial momentum, tangential mo- 
mentum, energy, and the equation of state, respectively, are 


dpi 

dr 


1 d 

+ ^("*0 


= 0 , 


(2.11a) 


dpi _ Pivl ZviVp 

dr f f 


~ dvp 

pi ir + 


dv i _*dv o u 0 a; 

+ + ~ 


d_ 

dr 


1 d_ 
r dr 


(™i) » 


(2.116) 


(2.11c) 
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9Ti ( 7 - 1 ) d 

+ r ^(ruj = 7 ( 7 - 1 ) 


dr 


dr 



x 7ia f.m 

-j — — r 

Pr r dr V dr 


( 2.11 d) 
( 2 . 12 ) 


where u\ — u\Re. Note that (2.12) implies that, for a slightly compressible vortex, 
the 0(M 2 ) density and temperature perturbations are of the same order of magnitude 
as the 0 ( 1 ) pressure variation. 

Using (2.11a, d), and (2.12) we can derive a partial differential equation for the 
evolution of T\ in terms of the known reference flow velocity Vq and pressure po 


dT 1 
dr 


— f (3) = ( 7 - i) te 

rPr dr \ dr J ^ ^ y dr 


vo V dp 0 

r ) dr 


(2.13) 


The equation set (2.11a-c), (2.12), and (2.13) will be used to to describe the evolution 
of the 0(M 2 ) quantities; the related boundary conditions are 


and 


II 

O 

v\ = 0 

fi = o, 

afi 

dr 

dr 

dpi 

= 0 , 


dr ’ 



f £&i — > 0, 


v\ 0 , 


r —> 00 


{ Pi 0, Ti 0, 


(2.14a) 


(2.146) 


( Pi->0. 

After the reference flow variables p 0 and u 0 are determined, one can solve for all of 
the 0(M 2 ) perturbation quantities. 



Chapter 3 

The Reference Flow 


One can completely determine the reference flow pressure and velocity from knowledge 
of the vorticity u) 0 . Therefore, we begin this chapter by discussing the reference flow 
vorticity equation. In the process of solving this equation we will find a group which 
leaves it invariant and identify a set of conserved quantities for the solutions. 


3.1 Derivation 


Taking the curl of (2.96) we obtain an equation for the diffusion of vorticity 

duo _ 1 d fduo\ 

dr f dr \ dr J ’ 


where 


Id 

w ° = 


(3.2) 


Here we assume that u 0 is bounded at f — 0, and that — » 0 as f -» oo. 

Based on the form of the integral invariants involved in self-similar solutions of the 
one-dimensional heat equation [5], we expect the following area-integrated, positive- 
integer moments of vorticity to be conserved 



(3.3) 


23 
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Consider the invariance of (3.1) and (3.3) under dilations in oj 0 , r, and r, and 
translation in r 

f' = : e a r, t' = e b r 4- (e b - l)r 2 , = e c £ 0 . 


The invariance of (3.1) and (3.3) requires that 6 = 2a and c = -(j 4 - 2) a. The group 
is 




e a r. 


r = 


e 2a r 


+ (e 2a - l)7j, 






(3.4) 


and a = 0 is the identity element of the transformation. The associated characteristic 
equations are 

dr dr d£n 

(3.5) 


r 2 (t 4- Tj) (j 4- 2)d>o 
Note that this group has three parameters: a, j , and r 2 . The resulting similarity 
variables are 

r 2 

71 = FT — V’ (3-6) 

4(r 4 - Tj) v 1 

and 

+ (3.7) 


Substituting these group invariants into (3.1), we obtain an ordinary differential equa- 
tion in the variable 77 

d 2 fl dfi (j + 2 ) _ . 

V &rf + ( V + ^drj + — 2 — ^ = °’ (° ^ V < °°)- (3.8) 

We require that Q is bounded at 77 = 0, and that H ■->> 0 as 77 -4 00 . Bateman [ 4 ] 
shows that the two-dimensional diffusion equation is satisfied by a separated solution 
in the variables r and 77 , and his solution method also leads to (3.8). For arbitrary 
values of j, he gives the general solution of ( 3 . 8 ) as 


= Cj exp(- 7 ?)M(-|; 1 ; t ;), (3.9) 

where Cj is a constant and M(— 1; 77 ) is a confluent hyp ergeome trie function. Note 
that an overbar and an upright font are used to distinguish this function from the 
previously defined Mach number M. 
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Using (3.6) and (3.7) the integral (3.3) becomes 


— t • C oo 

Aj = 4 (J+2 ^ 2 7T ^ ?7 i/2 ^(7?)dr/, (3.10) 


where i is a positive integer. Using (3.9) in this integral we find that Aj is unbounded 
for odd values of j? (see Appendix C). Therefore, we will only consider even values 
of j. Taking j - 2m, the confluent hypergeometric function reduces to the Laguerre 
polynomial 


M(— m; 1; rj) = L m (r)) = 


e v d m 
ml drj m 


m 

(re~”) = E 

k — 0 


m!(— rj) k 
(k\) 2 (m — k)\' 


(3.11) 


This solution for 0(r]) agrees with the vortex solutions found by de Neufville [13]. For 
reference, the first five Laguerre polynomials are 


1 1 

L 0 = l, Li = 1-7?, L 2 = - (?? 2 -4 j? + 2) , L 3 = — - (V - 9?? 2 + 18t? - 6) , 

U = Tij 4 - 16r) 3 + 72 n 2 - 96 r? + 24). 

In the even-integer case the conserved quantity Aj converges to a finite value for every 
m. Taking Aj -> A m for j = 2m, (3.3) reduces to 


f oo 

A m = 2ir ^ r 2m ti 0 rdr. (3.12) 

Each A m is the 2 m th area-moment of the m th vorticity solution. As pointed out by 
de Neufville [13], when m = 0 the total circulation of the fluid is invariant and when 
m — 1 the total angular momentum of the flow AA is invariant. These two invariants 
are a subset of the more general expression for the reference flow integral invariants 
(3.12) developed here. When m = 0, A m = A 0 , which is the total circulation of the 
reference flow, and for m = 1 the integral invariant is related to the total angular 
momentum A\ — —2 Ad. 
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In (3.9) let Cj C m for j = 2m. The constant C m is related to the integral 
invariant A m and is found by evaluating (3.10). Using (3.9) and (3.11), we can re- 
write (3.10) in the following form 


A 


4 (m+ i ) 7rC . 


ml 


r 

Jo 


oo r\m 

7] m 


d Tj m 


(i/V) dr?. 


(3.13) 


Integrating this equation by parts several times and then applying mathematical 
induction to the results, we can reduce A m to the following expression 

An = l o °° T, m e-M V . (3.14) 


The integral in (3.14) is the Gamma function [1], and for integer powers of m, C m 
and A m are related in the following way, 


C r . 


-A-m 

(-l) ro 4( m +i)(m!)7 r ‘ 


(3.15) 


Using(3.7), (3.9), and (3.11) the solution for the reference flow vorticity can be written 
as 


oo 0 = 


C m e V L m (rj) 

(t 4- Tj)( m+1 ) ’ 


(3.16) 


where C m /r/ m+1) corresponds to the amplitude (at f = 0) of the m th self-similar 
vorticity solution (when r = 0). The velocity and pressure corresponding to (3.16) 
are given in the following section. 


3.2 Physical description of reference flow solutions 

Each eigenfunction (solution for u> 0 ) satisfying (3.1) corresponds to a different radial 
profile of tangential velocity and pressure. Using (3.4) in (3.2) and (2.9a), we can 
find the transformation group for v 0 and p 0 


v' 0 = e - (2m+1)a 


ho, 


p’o = e- 2 < 2m+1 >“p 0 


(3.17) 
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These transformations leave the boundary conditions (2.10 a,b) invariant, and yield 
the following self-similar forms 

V(r)) =v 0 (t + Ti) (m+1/2 \ (3.18a) 

and 

~P(v) = Po(T + Ti) (2m+1) , (3.186) 

and the corresponding boundary conditions 


rj = 0 


< 


y = o, 


idV n 
n ^= 0 ’ 


77—^00 


V-+0, 

V->0. 


(3.19 a) 


(3.196) 


For each eigenvalue m, a completely self-similar solution can be found for the reference 
flow. The mathematical form of each solution may be used to describe the physical 
nature of the flow it represents. 


3.2.1 Tangential velocity 

Integrating (3.2), one obtains a solution for the reference flow velocity 


f 2C 0 

— [1 - exp(-r/)] , 


wo = S 


m = 0 


2C, 


[ ~f ? T + r .w eX P(- ? ?) PUnfo) - , m > 1. 


(3.20) 


When (m > 1), the term in square brackets in (3.20) is a polynomial which has one 
root at 77 = 0 and ( m — 1) roots for 77 > 0. Therefore, the tangential velocity for the 
m til eigenvalue will change direction (m — 1) times as one varies f from the origin 
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to f -* ex) (e.g. figure 3.16). Self-similar vortices corresponding to values of m > 2, 
would consist of (m — 1) concentric shear layers. 

3.2.2 Pressure 

From (2.9a) we see that dp^/dr > 0. Therefore, p 0 increases monotonically from its 
minimum value, at the center of the vortex (f = 0), to p 0 — 0, as f -» oo. Using 
(3.20) in (2.9a) gives p 0 


Po 


- r —dr = 

Jf r 


■Cl 


i 


°° ex P(~2>?) [L m {rj) - 


2(r + rj)( 2m+1 ) 

The pressure for the first few eigenvalues is 




dp. (3.21) 


Po = < 


cl 


1 + exp(- 2 ? 7 ) - 2 exp(-r 7 ) 


( r + r <) L 
C\ exp(— 2rj) 
4(r 4- Ti) z 


Cl exp (—27?) 


2 r) 


+ Ei (vj) - Ei(2?7) 



2 (r + r*) 7 

where Ei(?y) is the exponential integral 


7T 3n 5 
~8 + 16 


5 ^ II77 2 137 / 11 

72 36~ ~24 24~ 48 


m = 0 


m = 1 


m — 2 


m = 3, 


(3.22) 


/•oo p" -7 ? 

Ei (77) = / d?7. 

77 

The radial distributions of vorticity, tangential velocity, and pressure are plotted 
in figure 3.1(a-c) for the first three eigenvalues m, with C m = l,r = 0, and n — 1. 
Notice that Cjq = 1 at r = 0 for all three solutions shown. However, as m increases 
both the corresponding maximum tangential velocity, and maximum pressure drop 
decrease. 
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Figure 3.1: Radial distributions of (a) vorticity u 0 , (b) tangential velocity v 0f and (c) 
pressure p 0 , for C m = 1, r = 0, and = 1 , m = 0; , m = 1; , m = 2. 

3.2.3 Decay rates 

The solutions of the reference flow vorticity (3.16), tangential velocity (3.20), and 
pressure (3.21) can each be factored into two parts. The first part corresponds to an 
inverse power of r, and the second part is purely a function of 77. The 77 -dependent 
part of each solution governs the spreading of its radial profile by viscous diffusion, 
whereas the other factor governs the solution’s time rate of decay. 

For arbitrary values of m, the vorticity, tangential velocity, and pressure have the 
following inverse time dependence 

w 0 ocr" (m+1) , Do oc r _(m+1/2) , p 0 oc r~ (2m+1) . (3.23a, b, c) 

Note that the incompressible pressure decays like p 0 ~ Dg. 

In §2.2 the governing equations were scaled by k and v m , which are determined 
from the initial velocity profile in a vortex. However, the core size and maximum 
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tangential velocity in a vortex are functions of time. As a self-similar incompressible 
vortex decays, its core size l (radial location of maximum tangential velocity), grows 
at a rate governed by viscous diffusion: / oc ri Using this lengthscale and the 
time dependent maximum tangential velocity v m (r), we can show that the Reynolds 
number Re and the Mach number AT, at the core radius of a vortex vary in time 
as 

Re oc r~ m and M* cx T ~(m+i/ 2 )^ (3.24a, b) 

Interestingly, for m = 0 both the circulation (based on maximum tangential velocity 
Vm and core size l) and Re* are constant. If m > 0, both Re* and M* decrease in 
time — with the Mach number decreasing more rapidly than the Reynolds number. 
This time-dependent behavior is consistent with the underlying assumption that the 
Mach number of the vortex remains small. 

In the next section we will show that the inverse time dependence of the reference 
flow quantities has important implications for the form of the solutions we can use to 
model a self-similar vortex. 

3.2.4 Superposed solutions 

The Laguerre polynomials are a complete set of orthogonal eigenfunctions, and (3.1) 
is a linear equation. Therefore, an infinite series of Laguerre polynomials could be 
used to construct a description of any arbitrary vortex reference flow with uj 0 — > 0 as 
r However, as de Neufville [13] notes, the convergence of these polynomials 

is slow, and in order to approximate a particular vortex, a very large number of 
terms may be required. In the present work superposed solutions will be avoided, 
not because of their slow convergence properties, but because such solutions are not 
self-similar. In order to find similarity solutions to the 0(M 2 ) perturbation equations 
(2.11a-c), (2.12), and (2.13), the reference flow solutions (which form the forcing 
terms in these equations) must be self-similar. 

To see that superposed solutions are not self-similar, let a superscript enclosed in 
parentheses denote the m th vorticity solution and consider the sum of the ml th and 
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the m2 th vorticity solutions 

(ml) , (m2) C m iG ~ v l>mi{v) C m2 e v L m 2 (r}) 

0 0 (r + Ti)( mlJrl ) (r + r z -)( m2+1 ) 

Multiplying each term by (r + r^)( ml+1 ) we get 


(3.25) 


(r + r i )( ml+1 )(4 ml) +4 m2) ) = C rn ie~’ l L m i(r]) 

+ C m2 e-”h m2 (v) (r + (3.26) 

If the vorticity distribution in (3.26) is self-similar, all terms on the right hand side of 
the equation must be a function of rj alone. This is only true if ml = m2; therefore, 
superpositions of self-similar solutions are not self-similar. 

From a physical standpoint, a time-dependent quantity is self-similar if the spatial 
distributions of the quantity, at different instants in time, can be obtained from each 
other by the same group transformation. Each term (eigenfunction) in the solution of 
a vortex composed of superposed eigenfunctions will decay at its own, distinct, rate. 
The parts of the solution corresponding to large eigenvalues will become negligible 
long before the terms corresponding to small eigenvalues have decayed. 


3.3 The Oseen vortex, the Taylor vortex 
and Mandella’s vortex 

The similarity equations we will develop in chapter 4 are valid for any eigenvalue m. 
However, in §§5.1 and 5.2 we will focus on the solutions for the Oseen vortex (m = 0) 
and the Taylor vortex (m = 1). This section describes these two types of vortex, and 
qualitatively explores the formation of the vortex studied by Mandella [25]. 

The reference flow solution for m = 0 is the classical Oseen vortex [30], which 
is sometimes also called the Lamb vortex [20], or the Lamb-Oseen vortex [31], [15]. 
The solution when m = 1 is Taylor’s vortex [36] and has also been referred to as a 
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Taylor swirl” 1 [12]. Radial profiles of vorticity, tangential velocity, and pressure are 
given for both vortex solutions in figures 3.1 a-c. Note that, the tangential velocity 
profiles of both the Oseen vortex and the Taylor vortex have v 0 of the same sign 
(direction) at all values of f. However, as f -4 cx) the Oseen vortex has an infinite 
total angular momentum, and a finite circulation; the Taylor vortex has a finite total 
angular momentum and zero circulation. Also, note that the reference flow solutions 
for all higher order eigenvalues m > 2 have zero total angular momentum, and zero 
circulation for f — > oo [13]. 

The Oseen vortex is often used to model the viscous decay of vortex pairs such as 
those found far downstream from the wingtips of an airplane in steady, level flight. 
Assuming the vortices have the same circulation strength, but rotate in opposite 
directions, they move along parallel linear paths. Each member of such a vortex pair 
contains an infinite amount of angular momentum about its own center. However, 
the total angular momentum along a symmetry plane, equidistant from each vortex, 
and parallel to their motion, is zero — the circulation and angular momentum of each 
vortex in the pair are equal in magnitude, but opposite in sign. For the Oseen vortex 
Re is constant; see equation (3.24a). 

Now consider the vortex studied by Mandella [25], [26]. As shown in figures 3.2-3.3, 
this vortex is created by shock diffraction. A shock of strength M s ss 1.8 travels from 
left to right and passes over an airfoil. The airfoil has an angle of attack of 30°, and 
a thickness-to-chord ratio of 18%. How can the resulting flow be described? 

Currently there do not seem to be any reports, either computational or experi- 
mental, which give the vorticity, velocity, density, and pressure fields present in such 
a flow 2 . However, many features of the flow may be identified from the schemat- 
ics shown on the left side of figure 3.2. In these diagrams the heavy lines represent 
shocks, the thinner lines represent contact surfaces, the direction of propagation of 

X The use of this expression may have arisen to prevent confusion with the three dimensional 
vortices found in circular Taylor-Couette flow [14], Such vortices are also often called Taylor vortices. 

2 Moon &; Yee [29] use a TVD scheme to numerically determine the pressure and density of the flow 
for the early stages of Mandella & Bershader’s experiment [26]. Unfortunately, their computations 
do not capture the complete formation of the free vortex. 
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each shock is marked with arrows, and vorticity of positive and negative sign is de- 
noted by 4- and — , respectively. The accompanying holographic interferograms on 
the right hand side of figure 3.2 were provided courtesy of Dr. M. J. Mandella. 

The incoming shock first contacts the airfoil near the airfoil’s leading-edge. A 
reflected shock is produced which propagates away from the airfoil in a circular pattern 
(figure 3.2a). The incoming shock is split (into two parts) by the airfoil, forming an 
upper shock, which travels over the top surface of the airfoil, and a lower shock, which 
moves along the airfoil’s bottom surface. The circular, reflected shock is connected to 
the upper and lower shocks by Mach stem configurations. The upper shock diffracts 
and expands over the leading-edge of the airfoil and so becomes curved and weaker. As 
a result, the upper shock slows with respect to the lower shock, which is compressed 
and strengthened due to the airfoil’s high angle of attack. The lower shock reaches 
the sharp trailing-edge of the airfoil before the upper shock and diffracts around the 
tip. The diffraction produces a vortex of strong, predominately-positive (counter- 
clockwise) vorticity (figure 3.26). The curved upper shock and the small boundary 
layer on the upper surface of the airfoil both create negative (clockwise) vorticity. 
When the diffracted lower shock (now curving around to the upper surface of the 
airfoil) and the upper shock collide, the resulting shock structure continues to move 
from left to right. This shock structure then passes through the vortex at the airfoil’s 
trailing-edge, causing the vortex to shed from the airfoil and convect with the free- 
stream flow. In this free-vortex formation process, the strong, positive vorticity at the 
tail of the airfoil is encircled by negative vorticity of smaller magnitude (figure 3.2c). 

The schematic is highly simplified, and figure 3.3 shows that some elements of 
the real flow have been omitted: 1) The formation of a second, weaker vortex by 
shock diffraction at the leading-edge of the airfoil. This second vortex remains at- 
tached at the airfoil’s nose by a slip-line and seems to be accompanied by a small 
separation bubble. 2) A series of slip-lines or contact surfaces are present near the 
free vortex [26]. 3) When the upper and lower shocks collide near the trailing-edge of 
the airfoil, a focusing effect seems to occur. This focusing effect causes the resulting 
shock configuration to resemble the three-shock intersection structures reported by 
Sturtevant & Kulkarny [32]. This structure is visible in figure 3.3 as a triangular 
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(a) 



Figure 3.2: Three stages in the creation of a free two-dimensional vortex: (a) incident 
and reflected shocks, (b) diffraction around the airfoil’s trailing-edge, and (c) free 
convection of a vortex which contains both positive and negative vorticity. 
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Figure 3.3: Shadowgraph (left) and yellow & blue two-color schlieren (right) pictures 
showing the formation of Mandella’s vortex. The sequence is ordered from left to 
right and top to bottom. 
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system of contact surfaces and shocks (located slightly above and to the right of the 
free vortex). 

Despite their simplicity, the schematics in figure 3.2 capture one of the main fea- 
tures involved in the formation of a free vortex by shock diffraction: from physical 
reasoning, we expect such a vortex to contain vorticity of both positive and nega- 
tive sign. This idea motivates the examination of self-similar solutions for the Taylor 
vortex. Taylor’s vortex (m = 1), corresponds to the lowest and most stable eigen- 
function which satisfies the reference flow and possesses vorticity of both signs. For 
this reason, we will focus on analytical solutions for compressible perturbations to an 
incompressible Taylor vortex in §5.1. 

However, remember from § 1.1 that the flow in Mandella’s experiments is fairly 
complex, and that there may be acoustic waves present in the vicinity of the vortex. 
Again, the slightly-compressible vortex solutions developed in this thesis correspond 
to an idealized flow situation in which there are no acoustic waves and the far-field 
flow conditions are constant. Here, the Oseen vortex and the Taylor vortex are 
used as the base flows for the development of weakly- compressible vortex solutions. 
In a strict mathematical sense, there are shortcomings to representing a lone, two 
dimensional vortex with either the Oseen vortex or the Taylor vortex. If one assumes 
that the flow domain is completely unbounded and calculates the total circulation 
and the total angular momentum by performing integrations for 0 < r < oo, then: 
1) Oseen’s vortex has a non-zero total circulation but an infinite amount of total 
angular momentum — it would require an infinite amount of energy to create such 
a vortex, 2) Taylor’s vortex has a finite total angular momentum, but zero total 
circulation — such a vortex could not be produced by a lifting body, such as an 
airfoil, because it would necessarily shed vortices with nonzero total circulation 3 . In 

3 The lift on the airfoil and the distribution of circulation in the vortices in Mandella’s experiments 
has not been fully determined. As demonstrated by Taneda [35], the correlation between the lift 
on an impulsively started body and the presence of an isolated starting vortex is complex. He has 
shown that an impulsively started elliptical cylinder (2 : 1 aspect ratio), at angles of attack of 20° 
and 45°, experiences two peaks in lift magnitude. The first peak is associated with the development 
of a boundary layer and occurs in the absence of a visible starting vortex; the second peak arises 
just after the shedding of isolated vortices from both the leading- and trailing-edges of the body. 
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reality, these vortex solutions are only rough approximations of vortices existing in 
flows which can never be truly unbounded. 

These facts should be kept in mind when using the slightly-compressible vortex 
solutions developed in chapter 5. In the next chapter, the reference flow solutions are 
used to derive a set of equations which govern the self-similar forms of the compressible 
perturbations. 


Chapter 4 

Similarity Forms 

In §3.1 it is shown that the reference flow is invariant under the three-parameter 
group a, Ti , and j. Since the reference flow involves only even values of the integer 
j, j = 2m, we can alternatively use m as the third group parameter. It will now 
be shown that the compressible perturbation equations (2.11a-c), (2.12-2.13) and 
the corresponding boundary conditions (2.14a-6) are also invariant under this three- 
parameter group. We will use this group to find invariant forms of the 0(M 2 ) flow 
variables and the governing perturbation equations. 

4.1 Perturbation variables 

The similarity forms of the 0(M 2 ) compressible perturbation quantities are found 
by substituting the reference flow terms (3.7) and (3.18a-6) into the compressible 
equations (2.11a-c), and (2.12-2.13) and looking for a transformation group which 
leaves them invariant. We also require that this group yield the similarity variable 
77. Therefore, the 0(M 2 ) similarity solutions will have an f and r dependence of the 
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same form as the 0(1) similarity solutions. The relevant group is 


r = e a r, t' = e 2a r -f (e 2a — l)rj, $ == e _2(2m+1)a pi, ' 

Ul = e~( 4 m+ 3 ) a {t* ) p[ = e - 4 ( 2 m+l)a ^ i5 
v'y — e ~ 3 ( 2m + 1 ) a ^ 1 'p' — g~2(2m+l)aji^ 

Y = 7 , /V = Pr. 


(4.1) 


The gas properties 7 and Pr are treated as invariant quantities, and are assumed 
constant for this analysis. The characteristic equations of this group are 

dr _ dr _ dpi _ duj 

f 2 (r + Ti) 2{2m + l)pi ~ (4m + 3)u£ 


dgi _ dui _ dfi 

4(2 m + 1 )p 1 ~ 3(2 m + 1)77 ~ 2(2m + l)Ti ’ 


(4.2) 


and are used to find similarity variables corresponding to the 0(M 2 ) perturbed quan- 


tities 


<r(v) = ^(r + n)^ 1 ), 

P(w) = pi( r + TiY 2m+1 \ 
c* (77) = u[{t + Ti)( 2m+s / 2 \ > 
<t>{rj) = hi(r + rj) (3m+3/2 ^ 
i’i'n) = Pi(rP ri) (4m+2) . ' 


(4.3o-e) 
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4.1.1 Decay rates 

The decay rates of the 0(M 2 ) terms can be determined by inspection of (4.3 a- e) 

T\ r~ (2m+1) , pi oc r ~ (2m+1 ), u\ oc r -( 2 ^+ 3 / 2 ) ; (4.4a, b ? c ) 

OC r -( 3m + 3 /2) 5 pj oc r“ (4m+2) . (4.4d, e) 

Comparing (4.4a-e) with (3.23 a-c) we see that 2\ and pi decay at the same rate as 
the reference flow pressure (3.23c); this result could have also been inferred from the 
0(M 2 ) ideal gas law (2.12). The other 0(M 2 ) terms decay faster than the reference 
flow terms because the Mach number of the vortex decreases in time (3.246), and so 
the compressibility of the vortex also decreases. When m > 1, v x and pi decay more 
rapidly than the other 0(M 2 ) variables. 

4.2 Perturbed entropy, dissipation, and vorticity 

The solutions to (2.11a-c) and (2.12-2.13), form a minimum set of equations which 
can be used to completely describe the state of a compressible vortex — to 0(M 2 ). 
However, we can learn more about compressible vortices by explicitly investigating the 
entropy, dissipation, and vorticity within them. For this reason, we derive the equa- 
tions which govern these quantities before giving the invariant forms of the 0(M 2 ) 
perturbation equations (§4.3). 

4.2.1 Entropy and dissipation 

Using the Gibbs equation [37], the entropy (here in non-dimensional form) is related 
to the pressure and density in an ideal gas as follows: 



s = 


(4.5) 
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After writing p and p in non-dimensional form using (2.7), expanding (4.5) in a 
Rayleigh-Janzen expansion (2.8), and then approximating logarithmic terms as 

ln(x) = (x - 1) + i(x - l) 2 + • • • , (0 < x < 2), (4.6) 

[1] one can show [12]: 

s = (po~ Pi)M 2 + 0(M 4 ). (4.7) 

Hence, the 0(1) entropy variation of the reference flow is s 0 = 0 — which is expected 
since T 0 and p 0 are both constant. The 0(M 2 ) compressible entropy perturbation is 

h=Po-pi- (4.8) 


To find the dissipation, we take the total derivative of s, which is, to 0(M 2 ): 

Dsi _ ds 1 6 . ds 0 

Dt ~ dr +ReUl dr^' (4.9) 

=o 


Therefore, differentiating (4.5) by r gives the total derivative of the entropy, to 0(M 2 ), 
as 

Dsi _ dh _ dp 0 dpi 

Dt dr dr dr ' ( 41 °) 

Using (2.11a), (2.12), and (2.13), the local rate of change of entropy is given by [12] 


= (7 _ 1 ) (^s. _ hY + -LL (?A±) 

dr \dr f ) rPr dr y dr J 


(4.11) 


The first term on the right hand side of (4.11) is proportional to the rate of entropy 
production by viscous dissipation, and the second term is proportional to the rate of 
change of entropy due to heat conduction. In §6.4 we will use analytical solutions of 
the self-similar equations to investigate the relative magnitude of the terms in this 
equation. 

Before proceeding, note that (4.11) has an interesting consequence for the total 
heat flux in a slightly-compressible vortex. Using the fact that -» 0 in the far-field, 
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together with the 0(1) solutions for v 0 (3.20), one can show that the divergence of 
the heat flux into or out of the vortex in the far-field is zero, 


1 d 
Jim rxv 

r— >oo r or 



= 0 . 


Integrating this equation with respect to f gives 



where C is an arbitrary constant. In order for this relation to hold will vary as 


T\ ~ C ln(f) 


in the far-field. But, T\ — > 0 for f — > oo, therefore C = 0, and the far-field heat flux 
is zero, 


r ~ df 1 n 

lim r— — = 0. 

f— »■ oo or 


Because of this, the net heat flux within the vortex is also zero 




Regardless of whether or not a weakly compressible vortex is self-similar, if the overall 
temperature in the core of the vortex increases, the temperature rise is not caused by 
the conduction of heat from the far-field into the vortex. Instead, it is only caused by 
viscous dissipation and pressure work 1 . To state this in another way, the form of the 
perturbation equations limits the solutions to those cases for which one asymptotically 
obtains a uniform temperature (as the vortex decays) without the “feed” of heat from 
infinity. 


2 The total viscous work in a flow can be split into a part which is responsible for deforming fluid 
particles (the dissipation) and second part which accelerates fluid particles [37]. The acceleration 
part of the viscous work affects the kinetic energy in the flow, but does not increase the temperature 
of the fluid; see §6.2. 
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4.2.2 Vorticity 

Having determined the form of the 0(1) tangential velocity in §3.2 we can now derive 
an equation for the 0(M 2 ) vorticity. Expanding the tangential velocity in a Rayleigh- 
Janzen expansion (2.8), we get 


v — vq + M 2 V\ 0(M 4 ). 


(4.12) 


Taking the curl of (4.12) gives 


u — loq -f* M 2 ui -t- 0(./lT 4 ), 


(4.13) 


where 


Id 

^ = fd-f {rVl) 


(4.14) 


Upon taking of the curl of (2.11c) and also using (2.11a) we find the equation gov- 
erning the 0(M 2 ) vorticity 


dCb\ 1 d ( „du>i\ dpi _ ( dpi A du 0 _ du 0 

dr f dr \ dr ) dr W ° \ dr + Ul j dr Pl dr ’ 


(4.15) 


4.2.3 Similarity forms 

The similarity variables for the entropy, dissipation, and vorticity are 


CO?) = Si(r + T^ 2rn+1 \ 

Wv) = ^{^ + n) 2(m+1 \ > 

w(rj) = Coi(r + ri) (3m+2) , 


respectively. The corresponding 


Si a r 


— (2m+l) 


d§i 

dr 


decay rates are 

CC 7- -2 ( m+1 ) , Ul OC r “( 3m + 2 ) j 


(4.16a-c) 


(4.17a, 6, c) 
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and we see that the entropy variation decays at the same rate as T u p u and p 0 . 

In the following subsection the self-similar forms of each of the 0{M 2 ) equations, 
including (4.8), (4.11), and (4.15) will be derived. 


4.3 Perturbation equations 


Substituting the similarity forms of the perturbation variables (4.3 a—e) into equations 
(2.11a-c), (2.12), and (2.13), the partial differential equations describing the 0(M 2 ) 
perturbations are reduced to the following set of ordinary differential equations: 


d i d/3 . 

—( av *) = r 1 — + (2m + l)(3, 

= pV 2 + 2 4V, 


d 2 4> 


d(j) 




3 (m + -) — — 
/ 2 4r] 


(j)— (a ~ fir) 2 ) Q - m(3V, 
(3 = 'y'P -a, 


rj d 2 a (rjPr + 1) da 
Pr dr f Pr drj m ^ 


= (7-1) 


V 2 ( V' 
(2m 4- 1)V + — — ( C r 

1 \ rj2 / 


(4.18a) 

(4.186) 

(4.18c) 

(4.18d) 

(4.18e) 


Transforming (2.14a) and (2.146), gives the boundary conditions for (4.18a-e): 

i dip 


7 = 0 < 


a — 0 , 


0 = 0 , 


^ d^ = 0 ’ 


idcr i d/? 

,?2 ^ = 0 ’ ??2 4 = 0 ’ 


(4.19a) 


and 


?7 — > oo 


Oi -4 0, (j> — > 0, -0 -4 0, 

(7-4 0, (3 -4 0. 


(4.196) 
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The self-similar forms of the equations for the 0(M 2 ) density (4.8), dissipation 
(4.11), and vorticity (4.15) are 


C = -(7- l)V + a, 


(4.20 a) 


and 


<5 = ( 7 -i)(n-T) + 


J_d_ 
Pr drj 



(4.20 b) 


r) 


d 2 w 

dr/ 2 


+ (?? + !) 


dzu 

dr/ 


+ (3m + 2 )w 


(4.20c) 


The homogeneous parts of (4.18e) and (4.20c) have the same form as the similarity 
equation for the 0(1) vorticity (3.8). Therefore, their homogeneous solutions can be 
written immediately as 

a h = B m exp(-77 1 )L 2m (r? 1 ) (4.21) 

and 

^/i = W m exp(-77)L 3m+1 (?7), (4.22) 


where 1/1 = r/Pr , B m and W m are arbitrary constants, and the subscript h refers to the 
homogeneous part of each solution. Likewise, the homogeneous part of the equation 
for the 0(M 2 ) tangential velocity (2.11c) and the 0(1) tangential velocity equation 
(2.95) have the same form. Using (2.95), (3.18a), and (4.18c) one can show that the 
forms of 4>h and V(r + Tj)( m+ will be somewhat alike: 


<Ph 


(3m + l)Q m exp(—?7) 
1/ 2 


[L 3 m+l(»7) - L 3 miv)] , 


(4.23) 


where Q m is an arbitrary constant. 
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4.3.1 The behavior of the solutions 

Before giving solutions to (4.18a-e) for the Oseen vortex (m = 0) and the Taylor vor- 
tex ( Til 1) reference flows, we will use the forms of these equations to reveal some of 

the physical characteristics of slightly-compressible, self-similar vortices. Specifically, 
we will examine if these vortices can be isentropic, and briefly look at the Prandtl 
number dependence of the solutions. 

The similarity forms of the reference flow quantities constitute some of the forcing 
terms of (4.18a-e). The reference flow terms are functions of the viscous timescale 
r, and are independent of the Prandtl number. Pr appears explicitly only in (4.18e) 
— the equation for a. However, a and / 3 are directly related by (4.18<f), and either 
<T or S appear as part of the forcing terms in each equation of (4.18a— e). Therefore, 
all of the 0(M 2 ) terms will vary with both the viscous timescale r, and the Prandtl 
number as expected (§1). 

Can the compressible self-similar vortex be homentropic? In the homentropic 
case, si — C = 0. If this is true, then (4.20a) requires 

(7 ~ 1 = cr. (4.24) 

Using this relation to replace a with V in (4.18e), we see that V must depend on Pr in 
order for the vortex to be homentropic. However, since the base flow is incompressible, 
the base flow pressure must be independent of Pr. Therefore, a slightly-compressible, 
self-similar vortex cannot be homentropic 2 . The Pr dependence of the solutions is 
further explored in §§5.1.2 and 5.1.3. 

In the next section (4.18-4.23) are used to find self-similar solutions for the 0(M 2 ) 
compressible perturbations to both an incompressible Oseen vortex and an incom- 
pressible Taylor vortex. 


If a slightly-compressible vortex is not self-similar, one can specify an initial condition where 
the vortex is homentropic [12]. However, this condition is only valid at a single instant in time. 


Chapter 5 


Solutions 

5.1 Analytical solutions for the slightly- 
compressible Taylor vortex 

For m — 1 , (3.16), (3.20), and (3.22), can be used to reduce (4.18e) to the following 
relation 


d 2 a da 

dr)f dr/i 


H- 3cr 


= -Ci 2 (7~ l)exp(-27?) 


2 V 3 
r)~- + - 

2 4. 


(5.1) 


where, again, we take 771 = r}Pr. Analytical solutions to (5.1) for a general value of Pr 
can be found using the Method of Frobenius; the method is outlined in Babister [3]. 
The particular solution is 


a p = — C 2 (7 - 1)771 exp (—277! /Pr) 


D 0 + 2D 1 r }l + 4:D 2 r)l + £> n (277i) n 

71 = 3 


(5.2) 
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where 




A 


15 Pr 2 — 44 Pr -f- 32' 
144Pr 2 


(5.3) 


Ai = 


(4 - 2 Pr)D n -2 + [(2 n + 6)Pr 2 - 4Pr(2n + 1)] D n -i 
4(n 4- l) 2 Pr 2 


Using (4.21) we find that the homogeneous part of the solution for m — 1 is 

°h = — exp(-? ?1 ) (t? 2 - 4r)i + 2) , (5.4) 

and the complete solution is a = o h + a p . Using (3.22), (4.18J) and (4.20a), (5 and 
C are easily found. However, solving equations (4.18a-c) for a, and cj> is difficult 
because (5.2) contains a slowly converging series (see §5.1.2). 

Luckily, when Pr = 1, it is relatively easy to find simple closed-form solutions for 
most of the 0(M 2 ) terms. In this case rj = rj u and the particular solution to (5.1) is 

~ 7 4 ^ exp(— 2 v )(2t] + 1). (5.5) 

The total solution for cr is given by 

C 2 (y — 1) B 

° = exp(— 2 t7)(2tj + 1) + y exp(-j;) - 4 t? + 2) . (5.6) 

Using this result and (3.22) in (4.18d) and (4.20a) gives the solutions for and (: 

P = ex;p(— 2?y) [ 277(7 - 1) - 1] - ~ exp(-77)(?7 2 - 4 77 + 2), 


(5.7) 
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and 


C 2 h ~ 1) B 

C = — ^ 77 exp (—277) + ~ exp(—r))(r ] 2 — 4 t? -h 2). 


(5.8) 


Integrating (4.18a), we find the similarity parameter corresponding to the 0(M 2 ) 
radial velocity 


«= C?e ? ( i 2??) [(2~7) + 7/(1 - 2 7 ) + V( 7 - 1) 

4?72 


^772 exp(— 77) 


[77 s -677 + 6 


+ +■ (5.9) 

772 


Gi is a constant of integration which we will determine from the boundary conditions 
in (4.19a). 

The homogeneous solution for 0, which can be found using (4.23), is 


<t>h 


3- lV+36?? - 2 4), 


(5.10) 


where Qi is a constant. 

5.1.1 Initial/boundary conditions 

Ci and r< are determined so that £ 0 = 1, at r = 1, when r = 0: 

1 exp(l/2) 

* 2 ’ 1 2 

The radial velocity in the center of the vortex must be zero (4.19a), so 

G = g?(2-7) 


For this value of Gi, (5.9) shows that far from the vortex core, as 77 -+ 00, the radial 
velocity is 


Tii 


Gi_ C?(2- 7 ) 


77 : 


477: 


(5.11) 
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This result matches the prediction for the far-held radial velocity in the Taylor vortex 
made by Colonius et al. [12]. 

The solutions for cr, /?, and £ (5.6— 5.8) contain polynomials in rj. All powers of p 
in these polynomials are greater than zero. For this reason the boundary conditions 
(4.19a) are trivially satisfied for any value of B x . Here we will determine B 1 so that 
pi(f= 0, t = 0) is close to the value found in the experiments of Mandella [25]. The 
depth of the vortex density- well is defined as 


Ap* = 


p(r = 0,t)- Poo 

Poo 


— M 2 pi(r = 0, t). 


(5.12) 


If, at t - 0, the depth of the density- well is 55.0%, and the Mach number of the vortex 
is M — 0.67, then pi(f = 0 , r = 0) = -1.23. The forms of (5.2) and (5.5) differ, so 
in order for the series and exact solutions for the 0(M 2 ) density, temperature, and 
entropy perturbations to match (4.18d), (5.2), and (5.5) require that B Y for the series 
solution (denoted F?i s ), and B\ for the exact solution (denoted Bi e ) are related in the 
following way: 

B ls = B u - AllzF. (5.13) 

For an initial density-well depth of 55.0% and 7 = 1.4, B ls = -8.47 x 10 -2 and 
B\ e = -1.67 x 10~ 2 . In §5.1.2 we will show that the series solution for a (with 
Pr — 1) converges to the exact solution when a large number of terms are used 
in the series. The exact solutions for both the reference flow and the compressible 
perturbations are plotted in figure 5.1. 

The specification of the initial density-well depth is a departure from from the 
boundary conditions in (4.19a). Rather than only requiring that the gradient of the 
density is zero at f = 0, we are also specifying an initial value for the density at the 
origin. The boundary conditions in (4.19a) are invariant under the transformation 
group (4.1), but the initial density- well depth condition for pi is not. The resulting 
shape of the density distribution is dependent on the initial value chosen for Ap*. 
Since the coefficient Bi is determined by this initial Ap* the distributions of the 
self-similar solutions are dependent on the value of Bi and can differ substantially as 
Bi is varied. When B 1 is increased the temperature fi, at r — 0 increases and can 













CHAPTER 5. SOLUTIONS 


55 



Figure 5.1: Exact solutions for a slightly-compressible Taylor vortex with — 1/2, 
Pr = 1, 7 = 1.4, C\ = exp(l/2)/2, B\ e = —1.67 x 10 -2 . Results for four differ- 
ent times are shown: , r = 0; , r = 1/4; , r = 1/2, and , 

t = 1 . Reference flow quantities are plotted in (a) vorticity, (b) tangential veloc- 
ity, and (c) pressure. The compressible perturbations are shown in (d) temperature, 
(e) density, (f) entropy, and (g) radial velocity. 

become positive; the corresponding value of density pi decreases. The trends in pi 
and Ti are reversed when Bi is reduced (see figure 5.2). Since (4.19a) is still satisfied 
after applying the initial density-well depth condition, the solutions are self-similar 
in time. 


5.1.2 Convergence rates 

To determine the number of terms in the series solution IV, required to accurately 
give the distribution of the 0(M 2 ) quantities, 


T\ e ~ Tis 



err — 


(5.14) 
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is plotted in figure 5.3 for Pr — 1 and different values of 77 . The term err is a measure 
of the relative error in the series solution; T le and T\ s are the exact and series solutions 
for Ti, respectively. The error depends on both the size of N and the value of 77 at 
which the solution is evaluated. For larger 77 , more terms are required to accurately 
represent the solution (figure 5 . 4 ). 



Figure 5.3: Convergence of the series solution as a function of TV. Plots show the 
error (5.14) between the exact and the series solution for Pr — 1 , 7 = 1.4, t* = 1/2, 
T = 0, C x = exp(l/2)/2, Big = -8.47 x 10“ 2 , and B le = -1.67 x 10~ 2 . The error is 
calculated at three different locations: 

, 77 = 10 , , 77 = 20 , and , 77 = 30. 

5.1.3 Prandtl number dependence 

Consider how the 0(M 2 ) solutions vary with the Prandtl number (figure 5.5). Pr 
represents the ratio of viscous to thermal diffusion. When Pr — 1 , heat and viscosity 
diffuse at the same rate; when Pr < 1 , heat diffuses faster. Therefore, when the 
temperature profiles for different values of Pr are plotted on the same graph, and 
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Figure 5.4: Convergence of the series solution as a function of 77 for different N. Plots 
show the error (5.14) between the exact and the series solution for Pr = 1, 7 = 1.4, 
T i = 1/2, r = 0, Ci = exp(l/2)/2, B ls = -8.47 x 1CT 2 , and B le = -1.67 x 10~ 2 . The 
error is calculated for: , N = 15, , N = 30, and , N = 45. 

each profile is set to the same temperature at f = 0 and t — 0, the temperature 
distribution for the curve with the smaller value of Pr will extend further from the 
origin (figure 5.5u). Also, as Pr is lowered, the entropy distribution changes because 
the heat conduction term in the dissipation equation (4.11) varies inversely with Pr. 
As Pr decreases, heat conduction has a greater effect on the entropy variation s. 

In this thesis, values of Pr = 0.5, 0.72, and 1.0 are used to get a feel for the 
behavior of the 0(M 2 ) solutions when Pr is varied over an extreme range. Note that 
the Pr for most monatomic and diatomic gases at atmospheric pressure lies within 
the range 0.67 < Pr < 0.85 for temperatures between 100K < T < 1300K. For a 
given gas, the Pr is roughly constant, even at large temperatures. The theoretical 
value of Pr for a monatomic gas is Pr = 2/3 [37]. 
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Figure 5.5: Pr dependence of the series solution with r* = 1/2, r — 0, 7 = 1 . 4 , 
B\s = — 8.47xl0~ 2 , Ci = exp(l/2)/2. Radial distributions of 0(M 2 ) (a) temperature, 
(b) density, and (c) entropy, and are shown for three different values of Pr: 

, Pr = 1.00; , Pr = 0.72; and , Pr = 0.50. 


5.1.4 Entropy production 

In accordance with the Second Law of Thermodynamics, the total amount of entropy 
in the flow must increase. Far from the core of the vortex the entropy is constant, so 
for f — > 00 the entropy flux is zero. Therefore, the following inequality must hold for 
the total entropy in the vortex flow: 


d r°° 

— / sifdf > 0. 
dr Jo 


(5.15) 


As a partial check of the 0(M 2 ) perturbation solutions for the Taylor vortex, we use 
(3.6), (4.16) and (5.8) in integral (5.15) to verify that the Second Law is not violated. 
The result is 


d_ 

dr 



C?(7 ~ 1) 
2 (r + r ,) 3 ’ 


(5.16) 
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m 

c m 

n 

0 

0.699 

0.199 

1 

exp(l/2) 

2 

1/2 

2 

1.6402 

0.7964 

3 

4.5879 

1.0922 

4 

15.507 

1.387 


Table 5.1: Values of C m and n that give v 0 (r = l,r = 0) = 1 for the first five 
eigenvalues. 

which is always > 0 and so satisfies (5.15). Because the total far-field heat flux into 
or out of the vortex is zero (see §4.2.1), the net entropy produced in the flow is due 
to viscous dissipation. The magnitude of this dissipation depends on C Y . 

5.2 Numerical solutions 

The Oseen vortex (m = 0) and the Taylor vortex (m = 1) are both commonly used to 
model two-dimensional incompressible vortices. Therefore, it is interesting to compare 
the self-similar solutions of the compressible perturbations to these reference flows. 

The 0(M 2 ) perturbation equations (4.18a, d, e) and (4.20a) are easy to numeri- 
cally solve. As in §5.1.1 above, the reference flow constants C m and r* are fixed such 
that v 0 = 1 at f = 1 when r = 0. For reference, table 5.1 lists the values of C m and 
n which satisfy this condition for the first five eigenvalues. Equation (4.18e) is first 
numerically integrated to find the similarity parameter cr, which corresponds to the 
0(M 2 ) temperature variation. The second order Runge-Kutta Scheme is used to per- 
form the integration for 0 < rj < 10 on a 1000 point grid. The integration is started 
at rj = 0 and is initialized so that there is no the heat flux at the origin (dfi/dr — 0), 
and the temperature variation at r = 0 corresponds to a density variation with a 
density-well depth of 55.0%. The forcing terms in the equations are determined using 
the analytical solutions for the 0(1) velocity (3.20) and pressure (3.22). Equations 
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(4.3) and (4.16) together with (4.18a, d), and (4.20a) are then used to calculate u* v 
Pi, and si? respectively. 

The numerical solutions for an Oseen vortex reference flow with Pr = 0.5, 0.72, 
and 1.0 are given in figure 5.6. The numerical solutions of Ti, p u and s 1 for the 
Taylor vortex for these values of Pr duplicate the series solutions already shown in 
figure 5.5. Therefore, only the solution for u\ is shown in figure 5.7. 

The radial extent of both the Oseen and the Taylor vortex reference flows is larger 
for smaller values of Pr. However, as shown in figure 5.5 the temperature, density, 
and entropy distributions for the Taylor vortex exhibit peaks near the core radius 
at the smaller Prandtl numbers. The entropy distributions of both reference flows 
contain local maxima at the origin when Pr = 0.5 and 0.72, but entropy peaks at the 
core radius of the vortex do not appear when the reference flow is an Oseen vortex. 
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Figure 5.6: Pr dependence of the numerical solution for an Oseen reference flow, r* = 
0.199, r = 0, 7 = 1.4, Ci = 0.699. Radial distributions of C(M 2 ) (a) temperature, 
(b) density, (c) entropy, and (d) radial velocity are shown for three different values of 
Pr: , Pr = 1.00; , Pr = 0.72; , and Pr = 0.50. 
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Figure 5.7: Pr dependence of the solution for a Taylor reference flow. r £ = 1/2, r = 0, 
7 = 1.4, Ci — exp(l/2)/2. Radial distribution of the 0(M 2 ) radial velocity is shown 
for three different values of Pr: 

, Pr = 1.00; , Pr = 0.72; , and Pr — 0.50. 



Chapter 6 
Discussion 


One of the most significant differences between two-dimensional compressible- and 
incompressible- vortices, is that compressible vortices may possess a radial flow, but 
incompressible vortices cannot. Therefore, we begin our discussion of the 0(M 2 ) 
solutions given in §§ 5.1 and 5.2 by considering the radial velocity in a slightly- 
compressible vortex. 

6.1 Radial velocity 

The radial velocity u\ does not appear in either the 0(1) or the 0(M 2 ) radial mo- 
mentum equations, (2.9a) and (2.116), respectively. Therefore, the radial velocity 
does not arise from a force imbalance, such as a difference between the local pressure 
gradient and the centrifugal force in a rotating vortex flow. Instead, the form of 
the perturbation equations (2.11a-c),(2.12), and (2.13) and the 0(1) reference flow 
solutions (3.16), (3.20) and (3.21) requires the 0(Af 2 ) solutions to have a radial flow. 
The 0(1) reference flow pressure and velocity form the driving terms in the equation 
for the 0(M 2 ) temperature (2.13). The density p 1 is, in turn, found from the 0(M 2 ) 
ideal gas law (2.12), and so is dependent on T\ and po. The solutions for T\ and p$ 
necessitate the density in the vortex to increase with time. Therefore, as the vortex 
decays it is compressed, and the 0(M 2 ) continuity equation (2.11a) requires the ra- 
dial velocity to convect mass into the vortex. Starting with the 0(M 2 ) temperature 
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equation (2.13) we can show that far from the core of a slightly-compressible vortex 
“I is necessarily negative and varies like uj ~ 1/f. The volume integral of (2.13), 
taken over the whole flow-field, is 



Owing to the fact that there is no net heat flux into or out of the vortex (see §4.2.1), 
the second integral on the left hand side of this equation is zero. Using the 0(M 2 ) 
ideal gas law (2.12), (6.1) becomes 



The second integral on the right hand side of this equation can be expressed in terms 
of the 0(1) tangential velocity by integrating (2.9a): 


- t ~ [ Pordr = y rdr 

dr Jo dr Jo 2 

Using (6.3) and (2.11a) in (6.2) we find that 


(6.3) 


lim ru\ — 

f— too 




(6.4) 


Since the terms in the integral are positive definite, and physically realistic values of 7 
lie in the range 1 < 7 < 5/3, far from the vortex core the radial velocity is necessarily 
negative (see Colonius et al [12]). 

Note that as f — >■ 00 the compressible perturbation terms vanish (2.146), and so 
the far-field flow is essentially incompressible. Because of this, one would expect the 
asymptotic behavior of the radial velocity around the vortex to have the same form 
as the radial velocity around a point mass sink in a two dimensional, incompressible 
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flow. As we can see from (6.4), the far-field radial velocity has this property: it varies 
with f as u\ ~ 1/f. 

The far-field behavior of the radial velocity will be the same for all values of m 
and Pr , and suggests that a slightly-compressible vortex is always compressed. In 
figure 5.6d and figure 5.7 we see that for both the case of the Oseen and the Taylor 
vortex u\ is negative for all values of r. The magnitude of u\ is zero at f = 0, increases 
almost linearly to a maximum, and in the far-field decreases like 1 jf. However, just 
outside the vortex core (where r 1) the magnitude of the radial velocity decreases 
non-monotonically for Pr ^ 1. In both vortices, pi is larger near the core for the 
smaller values of Pr used, so the fluid within the core is compressed to a lesser extent 
— the corresponding radial velocity has a smaller magnitude. 

6.2 Vortex energetics 

The total amount of pressure work performed on a fluid volume is given by the integral 
of the product of the pressure and fluid velocity at the volume’s surface: 

J d ^ pu * dS = V • (pu) dV = J (u • Vp 4 - pV • uj dV. (6.5) 

Here we have used Stokes’ Theorem to convert the surface integral into a volume 
integral and split the total pressure work into two terms. The first term u • Vp 
represents the work required to move fluid against the pressure gradients in the volume 
and the second term pV • u is the work done by pressure forces in locally compressing 
the fluid at each point within the volume. In the far-field, u ~ 1/r, the pressure is 
constant, and the surface area of the two-dimensional flow field is proportional to r. 
Therefore, from the surface integral in (6.5) we see that the external pressure field 
surrounding a slightly compressible free vortex performs a net amount of work on the 
flow as the vortex decays. As viscous shear stresses decrease the tangential velocity, 
the centrifugal force of the flow’s rotation relaxes and the fluid in the vortex core is 
compressed by the external pressure field. This pressure work, in addition to heat 
released by the dissipation of the flow’s kinetic energy, increases both the temperature 
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and density of the vortex. To see this in more detail, we will consider the two parts of 
the total pressure work and their separate effects on the kinetic and internal energy 
of the flow. 

The scalar product of the velocity and the momentum equations gives an expres- 
sion for the balance of kinetic energy in a flow. Note that the kinetic energy associated 
with the radial component of the velocity does not enter our analysis because it is 
0(M 4 ). Integrating the product of v 0 and the tangential momentum equation (2.96) 
gives an expression for the total 0(1) kinetic energy in the vortex: 


d_ 

dr 





f dvp 

\ dr 



( 6 . 6 ) 


We see that since the reference flow is incompressible, only the work involved in 
viscous dissipation decreases the total 0(1) kinetic energy of the flow [21]. 

The product of vq and (2.11c), as well as the use of (2.9a), gives an equation for 
the 0(M 2 ) kinetic energy: 


Pi 


^(vj 

dr l 2 



+ ^ 


~*dpo 
- u 


i ^r +Vl df 


d .. _ , 


. d 


T d , ,1 

. d 

[Td 1 


+ f’O'XT 
dr 



(6.7) 


Integrating this expression over the entire flow field we find the total 0(M 2 ) kinetic 
energy in the flow 


d r°° 
dr Jo 



+ ViVq 


rdr 



1 df 


rdr 


\ or r ) 



rdf. 


( 6 . 8 ) 


From this equation, (3.6), (3.18a), and (4.3b, d), we see that the time derivative of 
the total 0(M 2 ) kinetic energy varies as l/r 4m+2 . Since the total 0(M 2 ) kinetic 
energy is a perturbation to the total 0(1) kinetic energy, it may be either positive or 
negative in sign. Since this term may be negative, the second term on the right hand 
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side of (6.8), which represents the contribution of viscous forces, is not necessarily 
positive-definite. However, in the case of both the Oseen and the Taylor vortex, the 
first term on the right hand side of (6.8), which gives the u • Vp part of the pressure 
work in (6.5), always increases the 0(M 2 ) kinetic energy of the vortex. The reason 
this happens is because the radial velocity in these two vortices is negative for f > 0. 

Turning now to the internal energy of the flow e, let e = C v T^e. With this choice 
of scaling, e = T and from (2. lid) we find an equation for the balance of internal 
energy in the flow 


dr \ dr r J Pr r dr \ dr J 


( 7 - 1 ) d 
r dr 


(ru[). 


(6.9) 


Integrating this expression over the entire flow field gives 


d_ 

dr 


r 00 /*c 

/ eifdf = 7(7 — 1) 

Jo Jo 


( dvo 

\ dr 



roo l Q 

ft - !) I 


( 6 . 10 ) 


note that the integral of the heat flux is zero (see §4.2.1). 

Because the first term on the right hand side of this equation is positive definite, 
and since fluid is moving into the vortex, the total internal energy of the vortex always 
increases in time. The second term on the right hand side of (6.10) is proportional to 
the 0(M 2 ) divergence of the flow and represents the total amount of work performed 
by the pressure in compressing the fluid within the vortex. To see this note that, to 
0(M 2 ), the term — pV • u in (6.5) is given by 


1 d 

_pv • u -4 (ful) , (6.11) 

where (2.7) has been used to make the terms non-dimensional, and the Rayleigh- 
Janzen expansion (2.8) was applied. 
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6.3 The entropy distribution 

For the initial and boundary conditions chosen in §5.1.1, the far- field entropy is larger 
than the entropy in the core of the vortex. What could produce a vortex with this 
feature? Once again, consider the formation of the vortex studied by Mandella & 
Bershader [26]. The experimental vortex isn’t necessarily self-similar, but here we 
will use Mandella s results to provide some insight into the entropy distributions 
in compressible vortices. Figures 3. 2-3. 3 show the formation of a vortex by shock 
diffraction. When the vortex forms, the shock passing beneath the airfoil diffracts 
around the trailing-edge of the airfoil (see §3.3) . The diffracting shock expands almost 
cylindrically around the sharp corner and becomes highly curved. During this process, 
the shock weakens substantially, and the entropy difference across it decreases. The 
vortex core forms at the trailing-edge of the airfoil, and is surrounded by fluid of higher 
entropy which wraps around the core as it rotates. The part of the shock passing 
over the upper section of the airfoil collides with the diffracted lower shock, and the 
resulting structure knocks the vortex away from the airfoil. The fluid surrounding the 
vortex has a larger amount of entropy, because the fluid comprising this region passes 
through the less-curved parts of the upper and lower shocks. Both the upper shock 
and the majority of the lower shock are stronger than the part of the shock which 
diffracts around the airfoil’s trailing-edge. Additionally, any portion of a shock which 
passes through the vortex, such as the shock structure which forms from the collision 
of the upper and lower shocks, is considerably weakened as it travels through the low 
pressure region of the vortex core. In this scenario one would expect the vortex to 
contain a local “cold spot” with less entropy than the surrounding fluid. If this is the 
case, does the experimental data reported by Mandella [25] support this conclusion? 

Mandella [25] presents formulae for curves fit to experimental measurements of 
pressure and density in a compressible vortex: 
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and 

A p 

P^Poo 7 7-2. (6.13) 

. . / r \ 


Here, Ar p = 0.117 cm is a length scale for the pressure profile, and A r = 0.166 cm is 
a length scale related to the density profile. The other parameters include: 

Poo = 1-787 kg m -3 , Ap = 1.076 kg m -3 , (6.14) 

and 

Poo = 152100 Nm" 2 , Ap = 107300 Nm~ 2 . (6.15) 

Also note that Mandella gives a core size « 0.93mm based on his estimate of the 
tangential velocity. 

We can use .these formulae to calculate the entropy s, in the vortex. Evaluating 
(4.5) with these parameters gives the curve shown in figure 6.1 d. For the values 
given in the above curve-fits, the entropy in the center of the vortex is larger than 
in the far-field. However, Mandella estimates that the pressure curve-fit is accurate 
to within ±5.6%. Note what happens when Ap is increased by 2.0% and then by 
4.0% from the value given in (6.15). The entropy in the center of the vortex decreases 
as Ap increases, and can become less than This demonstrates that the entropy 
distribution is very sensitive to small variations in the values of pressure and density. 
Therefore, one must be extremely careful to take into account the accuracy of the 
experiment. While the curve-fits may give reasonable measures of the pressure and 
density, they cannot be used to conclusively determine the exact values of s in the 
vortex; the prediction that the far-field entropy is larger than the core entropy lies 
within the accuracy of the experiment. 

6.4 Entropy production 

Most of the local entropy changes occurring in the core of a self-similar, slightly- 
compressible Taylor vortex are caused by heat conduction; since there is no net heat 
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Figure 6.1: Radial profiles of (a) pressure, (b) density, (c) temperature, and (d) en- 
tropy from curve-fits of experimental data (Mandella 1987). Small variations in the 
pressure are given by: 

, A p = 107300 N/m 2 nominal pressure curve; , A p = 109500 2.0% increase 

at f = 0, , A p — 111600 N/m 2 4.0% increase at f = 0. 
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flux into the vortex the total dissipation produced by heat conduction is zero. To see 
this, examine the dissipation equation (4.11). The rate of change of entropy is plotted 
in figure 6.2 for Pr — 1 at several different times. The maxima of ds\fdT correspond 
to the minima of §i (figure 5.1/). In figure 6.3 each of the terms in (4.11) is plotted 
separately. This allows us to compare the relative magnitudes of the rate of entropy 
changes due to viscous dissipation and due to heat conduction. In the core of the 
vortex r < 1, the magnitude of the heat conduction term is larger than the viscous 
dissipation term. Note that the part of ds\/&T due to heat conduction is larger than 
zero, because heat is flowing into the cold core of the vortex. For f« 1 and longer, 
the rate of entropy change from heat conduction is less than zero because heat is 
flowing away from this region and into the vortex core. The rate of dissipation from 
viscous diffusion is, of course, always greater than zero. As shown in figure 6.3, the 
series solution for <r p , (5.2), can be used with (4.16a) and (4.206) to reveal that the 
rate of entropy change from heat conduction increases in magnitude as the Prandtl 
number decreases; when fi is constant, heat diffuses faster, but the rate of viscous 
diffusion remains the same. 


6.5 Similarity and decay rates 

In §§4.1.1 and 4.2.3 the similarity forms of the 0(M 2 ) flow variables are used to 
determine the decay rates of the self-similar solutions. Barenblatt [5] shows several 
examples of physical phenomena which do not have self-similar initial conditions, but 
which asymptotically approach a self-similar solution at large times. Here we test 
this idea, and check to see if the self-similar solutions for the slightly-compressible 
vortex represent the limiting behavior of compressible vortices which are not initially 
self-similar. 

Consider the decay of the 0(M 2 ) flow quantities governed by (2.11a-c), (2.12), 
and (2.13). We will assume that the 0(1) flow is self-similar because, as discussed 
in §3.2.4, we cannot superpose the 0(1) reference flows for different eigenvalues m 
and still find self-similar solutions for the 0(M 2 ) terms. Since the reference flow 
is self-similar and the forcing terms in (2.11a-c), (2.12), and (2.13) are 0(1) flow 
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Figure 6.2: Entropy production in the Taylor vortex when Pr = 1, C m — exp(l/2)/2, 
n = 1/2. , r = 0; , r = 1/4, , r = 1/2, , r = 1. 

variables, the forced parts of the 0(M 2 ) solutions will also be self-similar. Note that 
these equations are linear in the 0(M 2 ) terms and the Laguerre polynomials, which 
constitute our solutions, are a mathematically complete set of functions. Therefore, 
after determining the corresponding particular solution, we can superpose the 0(M 2 ) 
homogeneous self-similar solutions to construct a solution for any 0(M 2 ) initial con- 
dition. At large times the vortex flows resulting from such a superposition process 
will be governed by the most slowly decaying self-similar component of the flow. 

For example, consider the decay rate of the 0{M 2 ) temperature fi from an initial 
condition which is not self-similar. Let the 0(1) base flow be a Taylor vortex m = 1, 
so that the self-similar part of the solution for the temperature T Us is given by (4.3a) 
and (5.6). We can alternatively represent this part of the solution by 

T - jlM 

^ (t + Ti ) 3 ’ 


(6.16) 
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Figure 6.3: Comparison of entropy production by heat conduction and viscous dis- 
sipation in the Taylor vortex when r = 0,C m = exp(l/2)/2,T* = 1/2: , viscous 

dissipation terms, and heat conduction terms for , Pr = 1, , Pr — 0.72, 

, Pr = 0.5. 

We add another self-similar term 


nh 

- 1 Ins — 


Mn) 

(t + n) ’ 


(6.17) 


to f lss so that the complete solution represents the evolution of the 0(M 2 ) temper- 
ature from an initial condition which is not self-similar 


Ti — T\ $ s + Ti ns = 


Mv) Mn) 

(t + Ti ) 3 (r + Ti) ' 


(6.18) 


This expression shows that the limiting behavior of the solution for fi is dominated 
by its most slowly decaying component, which in this case is Tins' 
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From (4.3 a-e) and (4.16a-c)we see that the decay rates of the 0(M 2 ) self-similar 
terms decrease as the eigenvalue m decreases. Therefore, the motion of a decaying 
vortex which is not initially self-similar, will asymptotically approach the self-similar 
solution corresponding to the component of the superposed solutions with the lowest 
eigenvalue. If a vortex under computational or experimental investigation is not 
initially self-similar, a self-similar solution will still give a close approximation to its 
motion at large times. 



Chapter 7 


Conclusions & Recommendations 


7.1 Conclusions 

The objective of this thesis is the development of a simple analytical description of a 
compressible free vortex. Applying group methods to the Navier-Stokes equations, it 
is found that they are not invariant under translations in the thermodynamic variables 
and so do not admit self-similar solutions for compressible free vortex flow. However, 
the equations governing a slightly-compressible free vortex possess this invariance 
property; these equations describe a set of incompressible reference flows and the 
compressible perturbations to each reference flow. The perturbation equations are 
derived under the assumptions that the viscosity and thermal conductivity of the 
fluid are constant, the far-field is quiescent with constant values of pressure, density 
and temperature, and that the flow evolves on a viscous timescale. The reference 
flow vorticity satisfies an eigenvalue problem in which each eigenvalue corresponds to 
a different self-similar incompressible vortex. An integral invariant for each vortex 
solution is determined. 

Using the perturbation equations it is shown that, generally: a self-similar slightly- 
compressible vortex cannot be homentropic; the far-field radial velocity in the vortex 
is negative and inversely proportional to the radial coordinate; and there is no heat 
flux into or out of the vortex in the far-field. 
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A three parameter group, which leaves the reference flow equations and the per- 
turbation equations invariant, is identified. This group is used to convert the flow 
variables to similarity form and to reduce the perturbation equations from PDE’s to 
ODE’s. The similarity forms of the flow variables are used to reveal their decay rates 
and it is found that all of the thermodynamic quantities decay like l/r 2m+1 . Addi- 
tionally, it is shown that superposed similarity solutions are not self-similar, and that 
slightly-compressible vortices will asymptotically decay to a self-similar state from 
arbitrary initial conditions. 

A self-similar solution for an Oseen vortex reference flow is found numerically, 
and both numerical and analytical solutions for a Taylor vortex reference flow are 
given. For general values of Pr the analytical solutions for the Taylor vortex are 
composed of infinite series. However, when Pr - 1 the solutions are simple closed- 
form expressions. The closed-form solutions are used to measure the convergence 
rate of the corresponding series solutions (with Pr = 1), and it is found that for 
0 < V IS 10 the sum of the first twenty- two terms of the series gives an approximation 
to the exact solution with an error of less than 1%. 

The solutions for the Oseen and Taylor vortices reveal that varying the Prandtl 
number has a significant effect on the distributions of temperature, density, entropy, 
and radial velocity in a slightly-compressible vortex. For Pr = 0.5, and 0.72 the 
radial temperature distribution for the Taylor vortex exhibits peaks near the core 
radius of the vortex. For these values of Pr, the magnitude of the radial velocity in 
both the Oseen and the Taylor vortex reference flows decreases non-monotonically 
from a maximum value near the core radius; in contrast, when Pr = 1 the decrease in 
the magnitude of the radial velocity is monotonic. It is shown that the non-monotonic 
decrease of the radial velocity’s magnitude, for Pr ^ 1, is due to local variations in 
the compressibility of these vortices. 

The density in the self-similar solutions is initialized so its minimum value in the 
vortex core corresponds to the minimum density in the vortex studied by Mandella 
[25] . With this choice of density at r = 0, the entropy in the core of the self-similar 
vortex is less than the entropy of the surrounding fluid. It is shown that this prediction 
for the entropy distribution lies within the accuracy of Mandella’s experiment and 
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that the entropy in the core of the vortex is extremely sensitive to small pressure and 
density variations. 

The dissipation in the Taylor vortex is split into two terms to show that most of 
the local entropy changes in the vortex core are caused by heat transfer. However, 
since there is no net heat flux into the vortex for r oo, heat conduction does not 
increase the total entropy in the vortex; instead it only redistributes the existing 
entropy. For the slightly-compressible vortex, the total entropy produced in the flow 
comes from viscous dissipation alone. 

An examination of the kinetic and internal energy in a slightly-compressible vortex 
shows that to 0(1) all of the kinetic energy in the vortex is consumed by viscous 
dissipation. However, to 0(M 2 ) the far-field pressure compresses the fluid within the 
viscously decaying vortex, thereby increasing the density in the vortex core. Part of 
this pressure work u • Vp, together with viscous dissipation, increases the internal 
energy and temperature in the vortex. The other part of the pressure work, pV • {?, 
increases the 0(M 2 ) perturbation of the vortex’s total kinetic energy. 

The search for self-similar solutions has yielded simple analytical expressions for 
the Taylor vortex, and a set of ODE’s which can be numerically integrated to find solu- 
tions for the compressible perturbations to a great number of different incompressible 
vortex reference flows. From a theoretical standpoint, the results of the above analysis 
provide a rich description of the basic physical properties of a slightly-compressible 
vortex and a foundation for understanding compressible vortices in general. From a 
practical standpoint, the use of the closed form Taylor vortex solutions is a quick and 
simple means of validating and initializing numerical simulations of flows containing 
a free two-dimensional compressible vortex. 

7.2 Future work 

A couple of important questions have arisen from this work and represent possible 
extensions of the current study. 

Here, it is assumed that the viscous diffusivity and thermal conductivity of the 
fluid are constant. In reality both of these quantities are temperature dependent. 
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An interesting continuation of the present work would be to include the temperature 
dependence of // and k by using power law models [41] 

fi — /i 0 T nl and k — n 0 T n2 ; 

for air n\ 2/3 and n2 — 0.81, and for Argon n\ — 0.72 and n2 — 0.73. Combined 
with a Rayleigh-Janzen expansion of the temperature, such models may permit one 
to derive a set of low Mach number perturbation equations which can be solved 
analytically. 

The analysis above does not account for the effects of acoustic waves in the flow. 
A theoretical study that couples the propagation of acoustic waves and the dynamics 
of a viscously decaying compressible vortex may help to explain some of the dis- 
crepancies between the findings of Mandella [25] and those of Colonius et al [12] 
(see Appendix A). Such a solution may also be useful for predicting/modelling the 
circular acoustic waves which arise in shock/vortex interaction [18]. 

Additionally, a full characterization of the flow during the formation and subse- 
quent evolution of the compressible vortex studied by Mandella [25], [26], has still not 
been done . Since compressible vortices are present in many flows of engineering and 
physical interest (see chapter 1), a thorough knowledge of the formation process and 
structure of a single compressible free vortex is of fundamental importance. There- 
fore, future computational and experimental effort should also be directed towards 
more fully understanding the flow studied by Mandella k Bershader. 


* 


1 See §3.3. 


Appendix A 

Angular Momentum 


In §6.1 it is shown that the radial velocity in a slightly-compressible vortex arises 
from the form of the perturbation equations (2.11a-c),(2.12), and (2.13) and the 
0(1) reference flow solutions (3.16), (3.20) and (3.21). Mandella & Bershader [26] 
suggest that the radial velocity in the vortex they investigated arises to counteract 
the outward spread of angular momentum by viscous diffusion. To further explore 
these ideas, consider the following preliminary investigation of the conservation of 
angular momentum in a compressible vortex. 

A dimensional analysis of the angular-momentum transport occurring in com- 
pressible vortices reveals three flow regimes into which their time evolution can be 
characterized. Each flow regime is determined by the magnitude of the radial veloc- 
ity and the core size of the vortex. Note that the main part of this thesis focuses 
on self-similar, slightly-compressible vortices, but the following analysis applies more 
broadly to compressible vortices in general. 


A.l Angular momentum transport 


For a two-dimensional and cylindrically symmetrical fluid flow, the equation governing 
the conservation of angular momentum, expressed in cylindrical coordinates, is: 


2 d fv\ d . . 1 

Pr Ft Ir ) + pU F (rV) ~ r 


ld_ 

dr 


' 3 d 
r ' W 

or 



(A.l) 
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The first term in this equation is the local time rate of change of angular momentum, 
the second its convection by radial mass flux, and the final term is its rate of change 
due to viscous torque. Using the constants li, r £ , p, u mj and v m , the angular 
momentum conservation equation can be cast in the following non-dimensional form: 



Pooh J 


i d r 3 a 

-WZ r 3 p~ 

r or dr 



= 0 . 


(A-2) 


For the analysis, numerical values characteristic of Mandella’s experimental results 
[25] are used for hi Poo, p, and v m (see table A.l). Note that v m cancels from the 
expression, and that the analysis yields the values for the time constant of angular 
momentum transport r £ , and the radial velocity (at h) u m . An examination of the 
orders of magnitude of the two resulting dimensionless groups (bracketed terms), sug- 
gests that there are three flow regimes into which the transport of angular momentum 
within a vortex can be categorized: convective, convective/diffusive, and diffusive. 


Constant 

Value 

Description 

h 

9.3 x 10“ 1 mm 

Vortex core radius 

Poo 

1.8 kgm“ 3 

Free-stream density 

p 

1.9 x 10“ 5 Nsm~ 2 

Coefficient of viscosity (air) 

V m 

2.3 x 10 2 ms -1 

Tangential velocity at 


Table A.l: Physical constants characteristic of the experimental vortex [26]. 


A. 1.1 Convective: T t u m /k ~ 1.0 > m/{Pooh 2 ) 

In this flow regime convection is the dominant process by which angular momen- 
tum is transferred within the vortex. By eliminating r £ in the inequality, approx- 
imate values may be found for the radial velocity: u m » 1.1 x 10 -2 ms -1 or 
more roughly u m > 3.0 x 10“ 2 ms -1 . Estimates for the time scale of this motion 
are given by r t ~ li/u m , which yields a time scale of r t <C 8.2 x 10“ 2 s or approx- 
imately r t < 3.0 x 10 -2 s. Note that the radial velocities predicted by Mandella 
u m ~ 0.3 ms 1 lie in the convective range, and that for the experimental vortex 
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under investigation, when u m is negative (inward flow) the angular momentum in the 
vortex is increasing. 

A.1.2 Convective/diffusive: r t u m /li ~ 1.0 - m/(Pook 2 ) 

When the magnitude of each dimensionless group is on of the order of 1.0, both 
the effects of viscosity and convection must be considered in determining angular 
momentum transport within the vortex. An approximate value of u m in this region 
is given by: 3.0 x 10- 3 ms' 1 < u m < 3.0 x 10" 2 ms" 1 . The radial velocities reported 
by Colonius et al. [12] u m ~ 6.8 x 10~ 3 ms -1 lie in this range. The time scale for 
these effects is r t ~ 7.9 x 10 2 s, which is the same as the viscous time scale of the 
incompressible reference flow solutions found in §3.1. 

A. 1.3 Diffusive: r t u m /li < 1.0 ~ (ir t / (pooh 2 ) 

In the diffusive regime viscosity is mostly responsible for the spreading of angular 
momentum. When slight compressibility is present radial velocities in the range 
u m < 1.1 x 1CT 2 ms -1 or around u m < 3.0 x 10~ 3 ms' 1 are expected for this regime. 

A. 2 Discussion 

The magnitude of the radial velocity measured by Mandella suggests that convection 
is the dominant mode of angular momentum transport in the vortex he describes. 
However, Colonius et al. suggest that the radial velocity should be around two or- 
ders of magnitude smaller than that determined by Mandella; in this case the angular 
velocity transport regime would be convective/diffusive. The discrepancy in the mag- 
nitude of the radial velocity may be due to unsteady variations in the experimental 
ambient flow conditions. In the computational study, the external pressure, density, 
and temperature are fixed — the only mechanism for changes to occur in the vortex 
is viscosity. In the experiment the vortex forms in a complex, unsteady flow field 
(§3.3). The vortex is created in the low-pressure region behind an airfoil, and after 
being traversed by a shock, is convected away from the airfoil. The radial pressure 
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gradient across the vortex is approximately balanced by its tangential velocity. If 
there are acoustic waves present in the flow, they could create small changes in the 
vortex far-held pressure. Such pressure variations might cause the tangential velocity 
profile to change, and the radial convection of angular momentum may be the mecha- 
nism responsible. Please see von Ellenrieder et al. [40] for a more complete discussion 
of how variations of the far-held how conditions may affect the transport of angular 
momentum and the distribution of entropy in a free vortex. 



Appendix B 

Introduction to One-Parameter 
Lie Groups 


This appendix gives a brief introduction to the basic ideas behind Lie (pronounced 
lee ) group methods. Much of the material presented here has befen adapted from 
Prof. Brian Cantwell’s course notes for AA218 “Introduction to Symmetry Analy- 
sis” [10]. A more complete description of the theory is presented in the course notes, 
and also in the text by Bluman & Kumei [8]. 

B.l Basic definitions 

A Lie group is simply a collection of parametrically-defined transformations on the 
variables in an equation 1 . An equation is invariant under a Lie group if the equation 
is symmetric under the set of transformations in the group. In geometry and art 
we say that an object is symmetric if we can perform an operation on it, such as a 
rotation or translation, and after the operation the object has the same appearance. 
In an analogous way, an equation is symmetric if we can substitute transformed 
variables into the equation and the resulting expression appears the same as the 

1 In general a group may also be applied to a system of equations. 
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original equation. For example, take the equation 

_ y 

dx x ’ 

and the transformation group 


(B.l) 


x = e s x, y = e s y, (B.2) 

where s is known as the group parameter. 

The differentials of x and y in (B.2) transform as 


d# = e s d:r, dy — e s d?/. 

On substituting the transformed variables into (B.l) we find 


e s d y e 5 y 

e~ s dx e~ s x 5 

and after canceling the exponential terms we get 


(B-3) 


d y _y 

df x 


(B.4) 


Aside from the tildes, (B.4) and (B.l) look the same; they are functionally equivalent. 
Therefore, (B.l) is invariant under the transformation group (B.2). 

A Lie group has four main traits, the group: 1) has an identity, 2) has an inverse, 
3) is invariant under composition (also known as “closure”) and 4) is associative. 
These properties are important in the development of self-similar solutions because 
they allow us to derive an invariant functional relationship between the variables in 
an equation. This functional relationship can then be used to simplify the equation. 

Let’s dissect the preceding paragraph by first examining the first four properties 
of a Lie group more closely. Returning to the example transformation group (B.2), 
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we see that when s — 0, e s = 1. So for a given value of the group parameter 5 , the 
group becomes an identity transformation. 2 

Next, consider a second transformation group for the transformed variables x and 
y in (B.2) above: 

x = e% § = e l y. (B.5) 

Rewriting (B.5) in terms of the variables x and y we get 


x = e* (e s :z) 
V = e t (e s y) 


= e t+s x , 
= e t+s y. 


(B.6) 


Here we see that if t = -5, then I = x and l = y. So the group (B.2) has an 
inverse which “takes you back” from the transformed variables x and y to the original 
variables x and y. Additionally, note that we can define a new variable u — t + s, 
such that 

~x = e u x, ~y = e u y. (B.7) 

The transformation (B.6), which is a composite of the transformations (B.2) and 
(B.5), appears the same as the original transformation (B.2) when expressed in terms 
of the ‘ dummy variable” u in (B.7). Thus, the group (B.2) is invariant under com- 
position. 

Now suppose that in forming the composite transformation (B.6) we had first 
transformed the transformation and then applied the result to x and y 


x = (eV) x — e t+s a;, 
V = (eV) y = e f+5 t/; 


(B.8) 


the composite transformation remains the same. Thus, the group is associative. 

To show how self-similar solutions to an equation can be found using a Lie group, 
let’s assume that an invariant functional relationship between the variables in an 

2 In general, a Lie group can always be mathematically manipulated so that when the group 
parameter is zero the group becomes an identity transformation. 
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equation exists. Using the vectorial representation of a Lie group 

x = X(x;s), (B.9) 

the invariant relationship between the quantities is given by 

#(x) = tt(X(x; s)) = T(x). (B.10) 

^ defines a family of characteristic surfaces or level curves such that T is constant 
when the invariant relationship between the variables is satisfied. The variables which 
compose the vector x then become parametric functions of s (figure B.l). 



Figure B.l: Along the invariant curve ip is constant and x and y are parametric 
functions of s. Lie group transformations map the point s 0 to the points s x and s 2 - 
£ is the vector of infinitesimals and is tangent to the curve ip. 

A convenient way of determining # is to arrange the group (B.9) so that s = 0 
corresponds to the identity transformation, and then expand the transformation in a 
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Taylor series near 5 = 0: 


s + 0(s 2 ). (B.ll) 

s=0 

The derivative in the second term, taken in the limit s -» 0, is known as an infinites- 
imal and is given by 


x = x + 


dX 

ds 


ax 


'dXdxi 


dx 

ds 

5 = 0 

dxi ds 

5=0 

ds 


(B..12) 


Here the subscript i denotes the i th component of the vector x. The sum of the first 
two terms in expansion (B.ll) is called an infinitesimal transformation 


X = x + fs. 

Substituting the infinitesimal transformation into (B.10) we get 

tf(x) = ^(x4-s£) 


= T(x) + s 


cM 

ds 


s 2 d 2 ^ 
T IT 


5=0 


2! ds 2 


+ 


d 3 V 


s = 0 


3! ds* 


+ 


5=0 


(B. 13) 


(B.14) 


With the relation 


(B.14) reduces to 


d 


'dXi dxi 

d 

ds 

i 5=0 

dxi ds 

s=0 dx ‘ 


(x) = (x) + s£ 


d$_ 

dx; 


- 2 d 


= ?«( x ) 


0$\ 


d 

dxi ’ 


2! ^ dxj y dx{ J 


5 3 ^ d 


+ & 


3! dx k 


c— (e — 

1 d'Xj V ' dxi 


+ 


(B.15) 


(B.16) 
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So a necessary and sufficient condition for <£ to be an invariant function is that 

(B.17) 


, dV „ 
65- = 0. 

OXi 


Geometrically, the gradient of each invariant surface in j s orthogonal to the corre- 
sponding vector field of the surface’s infinitesimals Since the gradient of a surface 
is parallel to the surface’s normal, the infinitesimals form a vector ( which is tangent 
to the surface (figure B.l). 

If there are n components in the vector of transformations (B.9), the partial 
differential equation (B.17) has an associated system of characteristic equations given 
by 

dx^ dxo dzs d.T^ 

(B.18) 


£1 £2 £3 £n 

The solutions of these equations determines the family of invariant surfaces T. There 
are n variables in the transformation group, and only (n — 1) characteristic rela- 
tions between them. Therefore, once the solutions ipj to the characteristic equa- 
tions are determined, a group of n variables can be specified by (n — 1) terms 
T = Wiipuifo , . . . , ipn-i)- As we will show below, this is where the power of Lie 
groups for simplifying and solving equations comes from. 

As an example let’s find # for the transformations given in (B.2). The infinitesi- 
mals are 


£1 = 


da; 
d s 


s = 0 


= x, and ( 2 -~ 
ds 


y- 


(B.19) 


5—0 


The associated characteristic equations are 


da; 

x 


V ' 


(B.20) 


and so the solution for T is given by 


T = 'L(^i), where ip x = 1 n(y/x) 


(B.21) 
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is a constant of integration. However ifj i, which is called an invariant , should not be 
strictly thought of as a constant, but as a parameter that determines which invariant 
curve the solution lies on. Along each curve the invariant will have a constant value. 
Let’s alternatively represent the family of invariant curves as 

T — 

If we let fa = y/x the relation between x and y everywhere on an invariant surface 
is given by 

V = fax. (B.22) 

Therefore, the family of surfaces 4/ corresponds to a set of rays centered at the origin 
(figure B.2). As we vary the invariant fa we move from one ray to another. But, 
along any given ray fa is fixed. 



Figure B.2: The transformation (B.2) maps points along rays centered at the origin. 
Each ray is specified by the numerical value of fa. 
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B . 2 Simplifying / solving equations 


In the simple example above there is only one dependent and one independent vari- 
able. Therefore, the relationship between x and y given by (B.22) is also the solution 3 
to (B.l). 

When there are more than two variables in an equation, the solution for T can be 
used to simplify the equation by reducing the number of variables required to specify 
a solution. In this case, the relation between the invariant surfaces can be rewritten 
as 

4>1 = f(lp 2, fa, • • • , fa- 1). (B.23) 

Take, for example, the linear, one dimensional, heat equation 


dcj) __ d 2 (f) 
dt dx 2 ' 


(B.24) 


This equation is invariant under the following Lie group 


<t> — fa x = e s x , t — e 2s t, 


(B.25) 


and the corresponding infinitesimals are 


fy = 0, fy = x, fy = 2 1, 


(B.26) 


respectively. Since </> is already an invariant, let fa = we don’t need to include it 
in the characteristic equations. The characteristic relations between x and t are 


drr dt 

x 2 1 


(B.27) 


Therefore, the second invariant relation is given by 


fa = xj \ft. 


(B.28) 


3 Note that when solving a problem in which boundary conditions are specified, the Lie group 
must leave both the equations and its boundary conditions invariant. 
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and the family of invariant surfaces is defined by 


^ = 'H'l’uih)- (B.29) 

As mentioned above, we can alternatively rewrite the relation for the family of invari- 
ant surfaces as Ui = /(^>)- Therefore, solutions to the heat equation can be found 
where 

4> = f(x/Vt). (B.30) 

If we set rj = x/Vt = ip 2 , the partial differential equation (B.24) can be recast as a 
very simple, readily solvable, ordinary differential equation 


d 2 (j> _ d (f) 

dr) 2 ^ dr) ' 


(B.31) 


One of the features of these solutions (B.30), which holds in general when the 
relation B.23 is used to simplify an equation, is that they are self-similar. In fact, 
one can generally say that the identification of a Lie group which leaves an equation 
invariant guarantees that the equation has a self-similar solution. 

Lastly, it should be noted that Lie groups provide an extremely powerful means 
of solving nonlinear differential equations. They are often used to identify integrat- 
ing factors for such equations. However, these topics are beyond the scope of this 
simple introduction; the interested reader should consult the references listed at the 
beginning of this appendix. 



Appendix C 

The Integral Invariants A? 


In §3.1 solutions for the reference flow similarity parameter Q(r]) and the correspond- 
ing integral invariants Aj are given for even values of the eigenvalue/stretching- 
parameter j. This appendix contains a more detailed explanation of the solution 
for 0 ( 77 ) and a proof that for positive, odd values of j, the integral invariants Aj are 
unbounded. 

Equation (3.8) is a form of Rummer’s equation [ 1 ] and its general solution is 


CM-^^ + DjUi-^v) 


Q{v) = e" 7 ? 

Here Cj and Dj are constants and U(—j/2; 1 ; 77 ) is 


(C.l) 


U(-^l-,ri) = A n 

2 sm( 7 r) 


M(— j/2; 1 ; rj) M(— j/ 2 ; 1 ; 77 ) 


r(-j/2)r(i) r(~i/2)r(i) J ’ 


which is mathematically “undefined”. Therefore, we take Dj = 0 and there is only 
one independent solution [3], which is given by (3.9). The hypergeometric function 
can be expanded to give 



H/2U n 

(1 ) n «! 


(C.2) 


97 


APPENDIX C. THE INTEGRAL INVARIANTS 


98 


( ) n is Pochhammer’s symbol [1] where 

(f) n ~ T(f) ~ + 2 ) • • • (/ + n — 1 ), 

/ represents any number, and 

(l)n = n! and (/) 0 = 1. 

Using (C.2) in the invariant integral (3.10) gives 

A . = 4 (H2)/2 wC , f tMn f°° ^/2 + n e - M)7 ( C . 3) 

n=0 Jo 

Evaluating this integral, we find that 

Aj = 4 (3+2)/2 <- E + " + !)■ (C.4) 

n=0 v 71 -; z 

First, let’s look at odd values of j. Both T(j/2 4- n + 1) > 0 and (n!) 2 > 0, but 
the sign of 

(~f )n = (~f )(-f + 1)(-| + 2) • • ■ + n - 1) 

depends on the value of n. For odd j, (~j/2) n changes sign only once as n increases 
from n = 0 to n — » oo. Because of this, the summation in (C.4) is unbounded. 

When j is even, (~j/2) n = 0 for n > j/2 + 1 and the series in (C.4) truncates. 
Taking j to be even gives a finite value of Aj and leads to the reference flow solutions 
found in §3.1. 



Appendix D 

Similarity at Higher Orders 


In chapter 4 the 0(M 2 ) compressible perturbation equations are shown to be invari- 
ant under the action of the same three-parameter Lie group (3.4) which leaves the 
reference flow equations invariant. In this appendix we will explore the similarity 
of the compressible perturbation equations for higher orders of approximation in the 
Rayleigh- Janzen expansion (2.8): 

/ = /o + M 2 fi + M A f 2 + M 6 fs + • • ♦ . (D.l) 


We will find that the 0(M 4 ) equations for the compressible perturbations are also 
invariant under the same group; the invariance of the perturbation equations may 
actually exist at each level of approximation in (D.l). 

To 0(M 4 ) the equations governing the compressible perturbations are 


dp2 1 d , . . . 

aT + faf (r “ 2 + ^“ lf) = 0 ’ 


(D.2a) 


dp 2 

df 


p 2 v 2 0 ( 2 piVQVi t 2v 0 V2 , v\ 

- j - I j — — 


~ dv 0 „ dm 

p2 ir +Pl ~dT 


dv 2 , (p^ + ul) d d , d 

+ 2 tS ("*>) + “ g= (™i) = Q-f 


dr 


dr 


1 d \ 
-fdf {rV2) 


(D.26) 


(D.2c) 
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df 2 . dfx 

F r + Pl 17 


+ <^- + (7-!) 


(pi +Ti) d 
f dr 


(ru{) + 


= 27(7 1) 


_- ]) fdv 0 vq \ f dv\ Vi 
\df r ) \df r 

IVi = P 2 + T 2 + pi7\. 


(7 ~ 1) d 


7 1 (9 
Pr f <9f 


Fr 



(D.2d) 

(M\ 

[ dr)’ 

(D.3) 


Here, each of the 0(M 4 ) terms is denoted by the subscript 2. 

The similarity forms of the 0(M 2 ) flow quantities are given in table D.l. Note 
that for each similarity variable, the next higher order similarity variable differs by 
a factor of (r + r .)( 2 ™+i). We will not pursue this pattern any further here, but it is 
quite possible that it persists at successively higher orders of approximation. 


0(1) 

0(M 2 ) 

0(M 4 ) 

V = Vo (r -h T*)( m ' f1 / 2 ) 
P = Po{ r + Tj)( 2m+ b 

a = f 1 (T + Ti )( 2rn +V 

P = Pi(r + nY 2m+1 ) 
a = u\(r + r-)( 2m+3 / 2 ) 

<j> = V\{r + ) (3m+3/2) 

0 =pi(r + ri) (4m+2) 

cr 2 — T 2 (r + Tj) (4m+2) 
02 = p 2 (r + ri)( 4m+2 ) 

a 2 = {^(r + Ti)( 4m+ 5/2) 
02 = 6 2 (t + rf)^ 5m+5 / 2 l 
02 = P 2 (r + n)^ 6m+3 l 


Table D.l: Similarity forms. 


The similarity forms of the 0(M 4 ) equations are: 

+ (4m + 2)0 2 ~ ^ ( a 2 7i l/ 2 + a/V /2 ) = 0, (D.4a) 


27? 


dip 2 

d?7 


= fcV 2 + 20V<I> + 2V<j> 2 + <i> 2 , 


(D.46) 
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-02 


’ dV ( 1 \J 

V d^ + { m+ 2) V 


- 0 


d(b ( 3 

T)-~ + f 3m + - 
dr] 


d(p 2 ( r 5 \ , 

V ^ + { 5m+ 2 ) ^ 


+ {0a + a 2 )^ (y /2 F) + a ^ (^V) 


drj 


= v 1/2 ^f {v 1 / 2 h ) , 


7 ?/j = 8 2 + a 2 + 0a, 


do 2 , , 


d a 

r) + (4m -f- 2) <72 

dr] 

- 0 

rj- — b (2m + 1 )cr 
[dr] J 


(D.4c) 

(D.4d) 


+ (7 - !)(<? + 0 )^ (»? 1/2 a) + (7 - 1) j- (j? 1/2 «2) 


27(7 - 1) 

' dV 

/n _ 

V' 

' d cj) 

4> 

V 

dr] 

" 2 

^ dr] 

2 


+ 2-1- (tip). 

Pr dr] \ dr] ) 

(D.4e) 

Note that even though each level of approximation in the compressible pertur- 
bation equations may be invariant, the sum of terms of different orders for a given 
variable will not be self-similar; this is because each level of approximation has a 
different decay rate. From a physical standpoint, we would not expect the sum of 
several different orders of compressible perturbations to be self-similar, because as 
discussed in § 2.1.1 above, the full Navier-Stokes equations do not admit similarity 
solutions for a free compressible vortex. 
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